LPNHE 2007-05 
FERMILAB-PUB-07-597-BSS 



The Dipion Mass Spectrum In e+e Annihilation and 
r Decay : A Dynamical (p, uj, cj)) Mixing Approach 

M. Benayoun^"), P. David^"), L. DelBuono^"), 
O. Leitner(«'^) and H. B. O'Connell^^) 

LPNHE Paris VPVII, IN2P3/CNRS, F-75252 Paris, France 
Laboratori Nazionali di Frascati, INFN, 1-00044 Frascati (Roma) Italy 
Fermilab, Batavia IL 60510, USA. 

March 18, 2008 

Abstract 

We readdress the problem of finding a simultaneous description of the pion form factor 
data in e+e^ annihilations and in r decays. For this purpose, we work in the framework 
of the Hidden Local Symmetry (HLS) Lagrangian and modify the vector meson mass term 
by including the pion and kaon loop contributions. This leads us to define the physical 
p, uj and (j) fields as linear combinations of their ideal partners, with coefficients being 
meromorphic functions of s, the square of the 4-momentum flowing into the vector me- 
son lines. This allows us to define a dynamical, i.e. s-dependent, vector meson mixing 
scheme. The model is overconstrained by extending the framework in order to include the 
description of all meson radiative {VP'^ and P77 couplings) and leptonic {Ve^e~ cou- 
plings) decays and also the isospin breaking (lo/c/) ^ vr+vr^) decay modes. The model 
provides a simultaneous, consistent and good description of the e'^e" and r dipion spectra. 
The expression for pion form factor in the latter case is derived from those in the former 
case by switching off the isospin breaking effects specific to e+e^ and switching on those 
for r decays. Besides, the model also provides a good account of all decay modes of the 
form VPj, -P77 as well as the isospin breaking decay modes. It leads us to propose new 
reference values for the — > e'^e~ and uj —>■ tt^tt^ partial widths which are part of our 
description of the pion form factor. Other topics (0 —>■ KK, the p meson mass and width 
parameters) are briefly discussed. As the e'^e" data are found perfectly consistent with 
r data up to identified isospin breaking effects, one finds no reason to cast any doubt on 
them and, therefore, on the theoretical estimate of the muon anomalous moment derived 
from them. Therefore, our work turns out to confirm the relevance of the reported 3.3 a 
discrepancy between this theoretical estimate of and its direct BNL measurement. 
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1 Introduction 



In order to study the phenomenology of light flavor mesons below 1 GeV, like partial de- 
cay widths or meson form factors, one needs a framework which includes in a well defined 
manner the lowest mass nonets of pseudocalar (P) and vector (V) mesons. Such a framework 
is well represented by the Hidden Local symmetry (HLS) model [T] [2l. In this approach, 
vector mesons are gauge bosons of a spontaneously broken hidden local symmetry which gen- 
erates their (Higgs-Kibble, HK) masses. Besides the non-anomalous sector, this model has an 
anomalous sector, hereafter called FKTUY Lagrangian [|3||, which aims at describing couplings 
of the form VVP, VP'j P77, VPPP for light flavor mesons. 

In its original form, the full (non-anomalous and anomalous) bare HLS Lagrangian fulfills 
the U (3) symmetry, as it possesses both Nonet Symmetry and SU{?>) flavor symmetry. As 
such, the HLS model covers a limited phenomenological scope. In order to broaden this scope, 
especially in order to account for VP'y and P77 couplings as derived from measured partial 
widths, the original HLS model should be supplemented with symmetry breaking mechanisms. 

Several SU (3) symmetry breaking schemes have been proposed [|4l El [H following the 
original idea of Bando, Kugo and Yamawaki (BKY) f7|. It has been shown [[810 that the most 
successful variant is the so-called "new scheme" of Ref . [[6l briefly recalled in Appendix O 
However, breaking only the SU (3) symmetry is insufficient flU |91 in order to reach a satis- 
factory description of the data on V P'y and P77 couplings which requires one to also break 
Nonet Symmetry. This was first performed in an ad hoc manner in [8] with the aim of recover- 
ing the radiative decay couplings of O'Donnell ifTOl which is the most general set fulfilling only 
SU (3) symmetry ; therefore, the model developed in (81 was indeed in agreement with gen- 
eral group theoretical considerations with additionally a mechanism smoothly breaking SU (3) 
flavor symmetry. 

Slightly later, it was shown ^ fTTl that an appropriate mechanism for Nonet Symmetry 
breaking can be produced by adding determinant terms [fT2| to the HLS Lagrangian ; the result 
was shown to meet all properties of Extended Chiral Perturbation Theory (EChPT) [fT3l[T4llT5l 
at leading order in the breaking parameters. Additionally, it was also proved [[TT| that the 
Nonet Symmetry breaking mechanism proposed in [81 was indeed an appropriate approxima- 
tion of the (rigorous) mechanism derived from adding the determinant terms to the bare HLS 
Lagrangian. 

However, breaking SU (3) and Nonet symmetries is still not enough to be in position of 
describing fully the whole set of radiative decays VP'y ; indeed, a process like (j) 7r°7 
requires including a — </> mixing scheme ; additionally, any global fit of all available VP'y 
transitions cannot be successful without introducing such a mixing. Traditionally, the to — (p 
mixing is described [8, 9] by rotating the fields luj and 0/ which are the entries of the bare 
vector field matrix, generally called ideal fields. An angle involving the mixing of the r] and 
7]' mesons is also required which has been shown IfTTI to vanish in the limit of exact SU{3) 
symmetry of the HLS Lagrangian. 

A brief account of the HLS model in its anomalous (FKTUY) and non-anomalous sec- 
tors is given in Appendices |A] and |Bl mostly focused on the subject of this work. Appendix 
O describes shortly but fully the various symmetry breaking procedures and the field renor- 
malization scheme, except for the SU (2) breaking mechanism which is the main subject of 
this paper. Appendix O deals with some specific features of the Yang-Mills piece of the HLS 
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Lagrangian which provides the kinetic energy part of the vector meson fields and should also 
undergo SU(3) symmetry breaking. 

These Appendices illustrate that the model depends on a very few number of free parame- 
ters as clear from Appendix |E] where they all appear. Among these free parameters, some are 
specific to our HLS based model, but some others are constrained. For instance, one of the two 
SU(3) breaking parameters is the ratio of decay constants fx / fn and can be constrained by the 
corresponding measured value. 

On another hand, decay processes like uj /(j) tt+tt^ can hardly be understood without 
some scheme for — uj and p° — mixings. Moreover, the decay behaviour of p° and 
cannot differ at the coupling constant level. It happens that the HLS model at one loop order 
provides a mechanism which allows us to perform a full — cu — (p mixing starting from 
the corresponding ideal (bare) fields. Within the non-anomalous HLS Lagrangian, the precise 
mechanism, kaon loop contributions, has been already described in [16J. However, we shall 
see that the (FKTUY) anomalous sector provides K*K loops as additional mechanism and that 
another one (K*K* loops) is provided by the Yang-Mills piece. 

Within the HLS model, as recalled below, Vi ^ Vj transitions among the ideal p% uji and 0/ 
fields are generated by loop effects ; these transitions generally have not a constant amplitude 
but rather depend on the squared 4-momentum flowing through the Vi (and Vj) line(s). We 
show in Section |4] that charged and neutral kaon, K*K and K*K* loops come through their 
sum in the uJi ^ (pi amplitude, while they come through their difference in p5 ^ 4'i and 
p5 ^ oji amplitudes. This also means that the uj — (p mixing proceeds from quantum effects 
rather than symmetry breaking effects, in contrast with the p° — or p° — mixings. These 
transition amplitudes are given by Dispersion Relations which should be subtracted in order 
to make the integral convergent (see, for instance. Appendix A in [9|). This gives rise to 
polynomials with real coefficients to be fixed using external renormalization conditions. 

In the exact SU(3) symmetry limit, charged and neutral kaons carry the same mass and one 
can expect (or require) the polynomials associated with the charged and neutral kaon loops 
to coincide. Likewise, the charged and neutral K*K loop functions can be made equal. In 
this case, the uji <-> 0/ amplitude survives with its renormalization polynomial, while the 
Pj ^ (pi and p5 ^ uJi transition amplitudes exactly vanish. If one breaks the SU(3) flavor 
symmetry leaving conserved the [u, d) quark sector, the same conclusion holds. However, if 
one introduces a breaking of SU (2) flavor symmetry, the mass difference beween u and d 
quarks generates a mass difference between the charged and neutral kaons (and K*'^). Then, 
the three possible transition amplitudes do not identically vanish any longer and, moreover, 
they depend on the invariant mass associated by the 4-momentum flowing through the vector 
meson lines. Stated otherwise : the p^, uji and cpi mixing into the physical p°, uj and cp fields 
should also be invariant mass dependent. As already noted in [|9lll6l this implies that the vector 
squared mass matrix, which has to be diagonalized in order to define the physical fields, is also 
invariant mass dependent and that the notion of mass for the physical vector fields becomes 
unclear as soon as one goes beyond tree level. 

The property that p^ — uj mixing should be invariant mass dependent has been the subject 
of several studies in the framework of general local effective field theories [17], then in Vector 
Dominance Models [[TSl [T9l where it was pointed out that the mixing amplitude should vanish 
at s = 0, as the p° and uj self-masses, in order to preserve gauge invariance. One may also 
quote other studies going in the same direction EOl |2T1| attributing the mixing to finite width 
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effects, or quark loops (and pion loops). Using only the pion form factor data available at that 
time, Ref. [ f22| derived an approximate expression for the p° — mixing amplitude ; however, 
limiting that much the kind of data used, one cannot really observe a clear mass dependence 
effect. On the other hand, one should also note that Ref. [1231 proved that isospin violation 
effects describing the — uj mixing are not accounted for by the low energy constants (LEC) of 
Chiral Perturbation Theory (ChPT) but are generated near threshold by the difference between 
charged and neutral kaon loops ; it was also shown that these effects are tiny at the two-pion 
threshold (10""^) while they are known to be at the percent level in resonance peak region. This 
illustrates the invariant mass dependence of the — uj mixing. However, one should note that 
the absence of the LECs in this calculation suggests a one-loop result may be unreliable. 

Therefore, the question we address is to consider the effects of loops on cu, mixing. 
However, in order to have some chance to single out the mass dependent behaviour, one clearly 
has to treat simultaneously the pion form factor data in the spacelike and timelike regions 
together with the largest possible set of light meson decays (radiative, leptonic and isospin 
violating strong decays) within a single framework. As argued above, the HLS model, equiped 
with suitable symmetry breaking mechanisms, seems able to provide such a framework. We 
shall not try to include Isospin Symmetry breaking effects into the coupling constants which 
would help by providing more parameter freedom in the fit procedures. If really needed, it can 
certainly be done along the lines of the BKY breaking scheme lITIO as illustrated by [241. 

The question of having a unified description of the largest possible set of low energy 
data is by itself interesting. However, this also addresses the puzzling and long-standing prob- 
lem of the difference between the (estimated) isospin 1 part of the pion form factor in e"'"e~ 
annihilations and in r decays which are related through the Conserved Vector Current assup- 
tion (CVC). The importance of the problem is enhanced by its implication for the predicted 
value of the muon anomalous magnetic moment to be compared with the direct BNL mea- 
surement [|25l ; referring to the latest account by M. Davier [,26 1 the estimate of the hadronic 
vacuum polarization (which enters the theoretical estimate) derived from e^e~ data provides a 
3.3 a disagreement between the theoretical estimation of and the BNL direct measurement 
[[25l ; moreover, the r data estimate of the hadronic vacuum polarization provides a value of 
very close to its direct measurement f26'[. 

Except for an experimental problem with e^e~ annihilation data (which seems by now 
unlikely) in the data recently collected at Novosibirsk [|27ll28ll29l[30l . or some new (or uniden- 
tified) physics effect, the disagreement between e+e^ and r data [l3Tll32l[33l is hard to explain. 
Indeed, a priori the single difference between these two channels, should be due to Isospin 
Symmetry breaking (IB) effects. However, the comparison has been already performed with 
IB effects accounted for in both e+e^ and r data. This includes [l34l [351 [361 pion mass values 
in kinematical factors, (a parametrization of the) p — uj mixing, charged and neutral p mass and 
width differences, short range W\ and long range [|38l [H |40l BB IH |43l IB effects in the r 
partial decay width to two pions. 

This persistent disagreement may point towards new physics effects ['441 ; however, one 
should also note that the way some IB effects are accounted for has been questioned several 
times. For instance, effects due to the charged and neutral p pole positions [[451 [461 were 
considered, but have not been found sufficient in order to solve the observed discrepancy [|35l 
[36l ; p — uj mixing effects may also have been poorly estimated ll47l . However, based on sum 
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rules derived using an OPE input, K. Maltman [|48l l49l concluded there is inconsistency of the 
(presently) estimated isospin 1 part of e+e~ data with expectation, while r data provide a nice 
agreement. 

We plan to address this question once again by building an effective model relying on the 
(symmetry broken) HLS model. In this approach, we plan to have a framework giving simul- 
taneously an account of the partial decay widths of light mesons decays (radiative, leptonic, 
isospin breaking decay modes), of the pion form factor in e+e^ data (both timelike and close 
spacelike regions) and in r decay. By construction, the corresponding expressions of the pion 
form factor will be such that they will solely differ from each other by Isospin Symmetry break- 
ing effects, mostly located in the (p, to, 0) mixing scheme which, of course, has no counterpart 
in r decay. Stated otherwise, our model is built in such a way that going from the pion form 
factor expression in e+e" annihilations to the expression valid for r decays is performed by 
switching off IB effects specific to e+e^ and switching on those specific to r decays. 

At start, the model we built is rendered complicated by the large number of possible loops 
involved. Fortunately, it can be somewhat simplified without loosing too much physics in- 
sight. Of course, this model depends on some parameters to be fixed in a fitting procedure ; 
we should define a fitting procedure flexible enough that r decay data can be removed or kept. 
Stated otherwise, the light mesons decays and the pion form factor in e+e data are expected 
to fix practically the U(3)/SU(3)/SU(2) breaking model. In this approach, we may get a predic- 
tion of the r decay 2-pion spectrum which can be compared with the existing measurements. 
Including r decay data should only refine the values of the fitted parameters. 

The paper is organized as follows : In Section [2l we derive the Lagrangian pieces of 
relevance in order to deal with the pion form factor in e"*"e~ annihilations and r decays, while 
in Section |3] we give the pion form factor expressions without loop corrections and symmetry 
breaking effects, mostly for illustration. In Section |4] we discuss the loop corrections which 
modify the vector meson mass matrix and perform already some simplifications. The modified 
vector meson squared mass matrix we propose is given in Section |5] with the diagonalization 
procedure and the relation between physical and ideal vector meson fields. The method used in 
order to renormalize the loop functions defining the self-energies and the transition amplitudes 
is sketched in Section [6l The form factor functions used for r decays and e+e^ annihilations 
are given in resp. Sections |7] and [U A necessary ingredient affecting the pion form factor in 
e"'"e~ annihilations is the photon vacuum polarization (VP), which is discussed in Section |9l 
Fitting with the partial width expressions is briefly discussed in Section [10]; more details can 
be found in [iT6l . where a practically identical method is used with constant mixing functions, 
however. The way to deal with the various kinds of data used in the fit procedures is described 
in SectiondU with a special emphasis on our dealing with some correlation phenomena present 
in the existing data. In Section [121 we fully describe the fit procedures we worked out under 
various conditions and the results and conclusions we reach ; we also comment on the numer- 
ical and physical properties of our model. Finally, Section [13] is devoted to a summary of our 
conclusions. 

As already commented upon in course of the Introduction, several Appendices gather re- 
sults of ours or others already published. They are given in an attempt to be as self contained 
as possible ; they are placed outside the main text in the interest of clarity and ease of read- 
ing. However, Appendix |D] contains new results concerning the Yang-Mills piece of the HLS 
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Lagrangian. 



2 The HLS Lagrangian Model 

We outline in the Appendices the main features of the HLS Model in both the non-anomalous 
(Appendix A) and anomalous (Appendix B) sectors. This allows us to derive the leading terms 
of the non-anomalous Lagrangian of interest for the present paper by expanding the exponen- 
tials defining the ^ fields. The breaking of flavour symmetries, SU(3) and the Nonet Symmetry 
is sketched in Appendix O 

Several pieces of the HLS Lagrangian of relevance for our problem will be given explicitly 
in the main text ; first, the part describing the photon sector (traditional VMD) is : 

vCvAfD = ^e(l - -n-dTr^ + i — izA-^- b)A ■ d K+ + i—hA -K^dlf 

2 Aza Aza 

(1) 

limiting oneself to vector mesons, pion and kaon fields. Flavor symmetries have been broken 
and, as noted in |l7l[6l, this implies a pseudoscalar field renormalization. The pseudoscalar field 
renormalization has been performed (following the prescription given by Eqs. (f86l) or, rather, 
by Eqs. (l90l) which include Nonet Symmetry breaking). The free Lagrangian of the vector 
meson fields is standard [1 , 2J, as well as the (canonical) pseudoscalar kinetic energy piece, the 
leptonic (see Eq.dTTl)) and photonic free Lagrangian pieces. 

The parameter g is the traditional universal vector meson coupling constant. On the other 
hand, the parameter a is specific of the HLS model and fulfills a = 2 in standard VMD ap- 
proaches ; however such a stringent condition is not mandatory and several phenomenological 
studies involving pion form factor data on the one hand [SOnSJJ and light meson decays on the 
other hand [81191 concluded that a much better favored value is a ~ 2.4 -^ 2.5. This opens a way 
to a direct coupling of photons to pseudoscalar pairs within VMD-like approches. One should 
remark the presence of a photon mass term of order which is traditionally removed by field 
redefinition [2l|52l (see also [53] and the discussion concerning the photon pole position). It 
can also be removed by renormalization conditions at one loop order. 

The parameter b in Eq. ([T]) is & = a{zv — l)/6 where zy is the SU(3) breaking parameter 
of the Cv part of the HLS Lagrangian, while za = [fK/Uf = 1-495 ± 0.031 (51 is the 
SU(3) breaking parameter of its Ca part HI IS. is almost fixed numerically, while zy is 
the major origin of the HK mass difference between the (j) meson and the (cu, p^) system ; it 
has to be fitted as the relation between vector meson masses determined experimentally and 
the (Higgs-Kibble) masses occuring in Lagrangians is unclear [[52l|9|. On the other hand, the 
value for fx / f-w provided by the Review of Particle Properties (RPP) ll54l can be added to the 
set of measurements to be fit. 



Pi + - 
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A subscript / on the fields, standing for "ideal", affects the neutral vector meson fields. It 
indicates that the corresponding fields occuring in the Lagrangian are not the physical fields. 

One should note |[6l that the SU(3) breaking of the HLS Lagrangian generates a non- 
resonant coupling of the photon to neutral kaon pairs ; this is a property common to all breaking 
procedures of the HLS Lagrangian proposed so far ||71|6l|4l|5l|. 

On the other hand, still limiting oneself to pions and kaon terms, the Lagrangian piece of 
relevance for r decay after symmetry breaking and field renormalization is given by : 



VudW^ ■ p- - 



lag 



a, 
2' 



za\Vi 



MS I ~ 



ud I 



7r° d TT" 



a 

+ 4I 



Za 



V2 



zv - Za\ 



IKsP 



(2) 

plus the conjugate of the interaction term (the term, not displayed). This Lagrangian piece 
depends on (72 (which is fixed by its relation with the Fermi constant (see Eq. (178]))). on the 
CKM matrix element = 0.97377 ± 0.00027 (Ml, on the universal coupling g and on the 
breaking parameters za and zy already defined. One should note, balancing the photon mass 
term in Cvmd-, a small mass term complementing the W mass of the Standard Model which 
could be removed by appropriate field redefinitions. 

Finally, the effective Lagrangian of the Model we use in order to describe low energy 
physics is : 

= ^VMD + Cr + Canomalous + C,Y M (3) 

where Canomaious IS givcn by Eq. (|931) . Its VVP part is not given in the Appendices, but can 
be found fully expanded in the Appendices of [8J. The last term is the Yang-Mills piece given 
by Eq. (|95] ) which undergoes SU(3) symmetry breaking as depicted in Appendix iDl: the main 
effect of this on the low energy phenomenology we deal with concerns the K* radiative decay 
widths (see the discussion around Eq. (|93])). 

The first two terms in Eq. ([3]) allow us to build up the pion form factor in e+e interactions 
and r decay and the leptonic widths of neutral vector mesons, while the anomalous decays will 
be dealt with starting from the third piece. All breaking parameters are common to all pieces of 
our C ; more precisely, all parameters of our model, except those of the vector meson mixing, 
could be fixed from only Canomaious and the leptonic decays of vector mesons. This was proved 
in |[8l|9lli6l by adding various (p, uj, (p) mixing schemes including the most traditional {cu, (p) 
mixing in isolation [[Si. The case of the parameter zt, which is important only for K* radiative 
decays, is special ; this comes out as the SU(3) breaking parameter of Cym which shows up in 
the anomalous Lagrangian of Eq. (|93l) when replacing there the bare vector field matrix by its 
renormalized partner (see Appendix iDl). 
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3 The Pion Form Factor Without Symmetry Breaking 



The Lagrangian given in Eq. ([3]) allows us to construct the pion form factor in e^e anni- 
hilation and in r decay. One has (m^r = m^±) : 



""^^ ^ 2' 2 Dv{s) 
for both processes involving intermediate photon and W boson. We also have 



(4) 



0" e ' e 



— (s) 
ds 



TT ' TT 



F^{s)\^[G,{s) + ^G2{s)] 



(5) 



wit 



itlfl: 



Go{s) 
G2{s) 



4 K - s)'K + 2g) 3 



,3/2 



(m^ - s)^(4s - m^) 4 



(6) 



,5/2 



and : 



mio — m'i+ 

7r TT^ 

2 



2 ; s - Ami 



-0.06455 



(7) 



Qn 



2v^ 



where one accounts for the pion mass difference. The (s) term gives a completely negligi- 
ble contribution to the form factor and will be cancelled out from now on. The bare inverse 
propagator Dy = s — has to be modified for self-mass effects which fortunately shift the p 
meson pole off the physical region by giving it an imaginary part. At this stage, there is also no 
inclusion of loop effects in 7p or Wp transition amplitudes in the expression for the pion form 
factor itself. 

Additionally, there is clearly no interplay of the to or (p mesons as can be seen from in- 
specting the various pieces of the full Lagrangian in Eq. ([3]) ; this should come from Isospin 
Symmetry breaking. 

Including self-mass effects for the p (and adding the co/cj) meson contributions for e+e^ 
annihilation), these expressions provide the usual HLS based framework for pion form factor 
fitting of e^e^ data fSTllSOll. Even if never done, in principle, Eq. & applies to r data, again 
after shifting the p"^ singularity off the physical region by means of a varying width Breit- 
Wigner amplitude, for instance. 



'Of course, in the SU(2) symmetry limit, we have e = and g^r = Q-k- 
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4 Including One-Loop Effects In The HLS Lagrangian 



From the expressions given in the previous Section, the p meson occurs as a pole on the 
physical region ; this is moved off the real axis by self-mass (loop) effects which essentially 
turn out to provide a width to the p through the imaginary part of the pion loop. However, 
besides this effect, a closer look at our C allows us to see that loop effects contribute to generate 
self-masses to all vector mesons, and transition amplitudes among all neutral vector mesons. 
Assuming from now on SU(3) and SU(2) breaking effects, the charged and neutral pion and 
kaon masses become different. One can see that kaon loops and the anomalous VVP piece of 
the Lagrangian give the following transition amplitudes for neutral vector mesons (ideal fields 
are understood, i.e. Yippis) should be understood as lip pis), for instance) : 



glJi'{s)+glKj,{Ii+{s) + Iio{s)) + 



n..(s) = ^?5x;,(n+(s) + no(s)) 



+ 



+ 



+ ■ 



+ 



91k*k^k*+k-{s) + ■ ■ 



^u4>{s) = -gwKKg4,KKiJ^+{.s) + no(s)) + [2 g4,K*KguiK'KiJ^K'±K^{.s) + n/^*oj^o(s)) 



gpKKgi^KKO^+{s) - no(s)) 



^p<i>{s) = -gpKKg4,KKiJ^+{s) - no(s)) 



+ [2 gu^K*KgpK*K{^K*±K^{s) " n^.o;^o(s)) 
+ [2 g<j>K*KgpK*K{^K*±KT{s) - n^.Oi^o(s)) 

(8) 

where we have defined gp^^ = ag/2, gpKK = g^KK = ag/ (Aza) and g^KK = V^agzy/ {4:Za). 
n'(s) is the charged pion loop, while n+(s) and no(s) are the charged and neutral kaon loops, 
each amputated from their couplings to vector mesons {i.e. loops carrying unit coupling con- 
stants). 

The contributions of the anomalous loops have been displayed between square brackets. 
The anomalous FKTUY Lagrangian gives several terms contributing to the self-masses Ilpp{s), 
^uiuj{s) and 1100 (s). They can easily be constructed from the VVP Lagrangian given in Ap- 
pendix 4 of flSll ; for these, we have displayed in Eqs. ([8]) only one representative of the full list 
which includes always K*'^K'^, K*'^TC^ ,K*^ and, depending on the particular self-mass 
considered, contributions from cuir^, p7r°, cur], pi], (pi], cui]', pi]' or (f)i]' loops. 

The anomalous parts of all transition amplitudes have been entirely displayed, as they ex- 
hibit an interesting correspondence with the non-anomalous contributions. We have identified 
to each other both K*'^K^ loops on the one hand, and K*^TC^ with K*^K^ on the other hand. 



Using the present set of notations, we have defined gpK-K = g^K-K = y zt/ ZAGanomf^ 
and g^K'K = Ganom/ V^zaZt with (HI Ganom = Sg^ / (Stt^ f^) . We have also denoted by 
IIk*±k^{s) and IIk*oko{s) the amputated K*^K^ and K*^K^ loop functions. 

In the exact Isospin Symmetry limit, one has n+(s) = no(s) ontheonehand, andnx*±i<'=F(s) 
n^.oxo(s) on the other hand. Then, all transition amplitudes vanish except for ujcj). 

Therefore, within the HLS model, the ujcf) mixing is a natural feature generated by loop 
effects and not from some breaking mechanism. In contrast, the p0 mixing and the prominent 
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puo mixing are fully due to Isospin Symmetry breaking. Including the anomalous sector does 
not change the picture. 

Actually, as emphasized in Appendix O the Yang-Mills sector of the HLS Lagrangian 
Il2l provides a further mechanism which comes in supplementing the kaon and K*K loop 
mechanism just described. This additional mechanism is produced by charged and neutral 
K*K loops ; these still come in by their sum in the ujcj) transition amplitude, while it is their 
difference which takes place in the puj and pcj) transition amplitudes. 

One may wonder that, within VMD-like approaches, all mechanisms contributing to the 
vector meson mixing at order always involve loops with a pair of mesons carrying open 
strangeness. This was true for the non-anomalous HLS Lagrangian (KK) and for the anoma- 
lous HLS Lagrangian (K*K) ; we also find it true for the Yang-Mills piece (K*K ). 



If one denotes by n+/o(s) the amputated K^K^ loop and by n"(s) the amputated vr 
loop, the charged p self-mass reads : 



±^0 



n;p(s) = 9l^Jl"{s) + 2gl^^Tl^„{s) + [gl^^^AAs) + ■■] (9) 

with a partial display of the anomalous loop contributions between the square brackets. The 
Yang-Mills term introduces also pp loops (not displayed as they are of little importance at our 
energies). This expression actually differs little from the neutral p self-mass ; indeed, the effect 
of having different masses for neutral and charged particles in these loop computations is tiny. 
Of course, in the Isospin Symmetry limit, we have n^^(s) = T\.pp{s) and then the p° and 
propagators (and their poles) coincide. 

As is clear from the Lagrangian pieces given by Eqs. ([T]) and Q, the fields p^, uj^ , (p^ 
as well as are certainly mass eigenstates at tree level. This statement remains true for p"^ 
at one-loop order as there is no transition loop from this meson to any other one. This is, 
however, clearly not true for p\ uj^ , cj)^ fields which undergo mixing with each other, as can be 
seen from Eqs. ([8]). Moreover, as for the self-masses, these transition amplitudes are invariant 
mass dependent as already noted [161 ' 

This implies that physical fields associated with the p°, u, mesons do not coincide with 
their ideal combinations as soon as one-loop corrections are considered. Moreover, the precise 
content of the physical fields in terms of ideal fields varies with s, or more precisely with the 
invariant mass flowing through the physical field under consideration. This does not prevent 
in standard approaches to use p^, uj^ , (j)^ in physical amplitudes Il55l . However, as one loop 
effects have certainly to be considered even only in order to shift the vector meson poles off 
the physical region, they should legitimately be considered also for field mixing. 

We raise the question of taking these loop effects properly into account and proceeding to 
the appropriate field redefinition to physical fields. In order to deal with this problem, let us 
define as effective Lagrangian the Lagrangian in Eq. ([3]) supplemented with the self-masses 
and transition terms occuring at one-loop order ; this turns out to replace the simple vector 
meson mass term in the HLS Lagrangian by (m? = ag'^f^): 

i {[m^ + Upp{s)]pj + [m^ + U^^{s)]uj + [zym^ + U^^{s)](Pj 

^mass \ (10) 

+2Ylp^{s)piUJi + 2Ilp^{s)pi(t)i + 2n,^(s)cjj0j} + K + n;,^(s)]p+p- 
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The K* mass term, which should also be modified correspondingly, is not shown as it plays 
no role in the present problem. 

Even if anomalous VP contributions seem to play some role visible [ISTTl (and, nevertheless, 
marginal) in pion form factor data, qualitatively their explicit form is really active only above 
the iun threshold, which is the lowest mass VP threshold ; all others are far above the GeV 
region^. Below the threshold, the main effect is due to their subtraction polynomials which can 
be well absorbed in the subtraction polynomials of the accompanying pion and kaon loops in 
order to put the poles of the p propagator at the place requested by the data. 

Beside the (non-anomalous) pion and kaon loops, all transition amplitudes involve K*K 
and K*K loops, the thresholds of which being resp. at ~ 1.4 GeV and ~ 1.8 GeV. This means 
that, besides their subtraction polynomials (minimally of degree 2), in the region below the 
GeV, their behavior flU is a real logarithmic function (below ^/s ~ 0.4 GeV) or an arctangent 
function (0.4 < y/s < 1.4 GeV). This also can be numerically absorbed in a fitted subtraction 
polynomial. 

Therefore, there is some sense in neglecting the (explicit) contributions of the anomalous 
and Yang-Mills loops, being understood that their effect is mostly concentrated in their sub- 
traction polynomials. Moreover, as these come always together with pion and kaon loops, they 
can be accounted for by simply letting the (free) pion and kaon loop subtraction polynomials 
to be second degree. Therefore, we shall neglect their (explicit) contributions, pointing at the 
appropriate places to their possible influence. Thus, from now on, the self-energies and transi- 
tion amplitudes should be understood as Eqs. dHHH) amputated from the square bracket terms 
and without the Yang-Mills contributions depicted in Appendix iDl 

The use of the modified HLS Lagrangian has already been discussed in (9) where it was 
shown, for instance, that this method allows one to recover vector meson propagators usually 
derived through the Schwinger-Dyson resummation procedure, which turns out to sum up an 
infinite series, which is not necessarily convergent. However, we show shortly that introducing 
this modified mass term allows us to also account for the other transition effects which would 
be more difficult to derive from the Schwinger-Dyson resummation procedure (of course, this 
should be possible, merely tedious). 

5 Mass Matrix Diagonalization And Physical Fields 

As clear from Eq. (fTOl) . at one loop order, the mass term is diagonal in the charged vector 
meson sector and will not be discussed any longer. In the neutral vector meson sector, however, 
the mass matrix is not diagonal and the effective Lagrangian mass term is : 

Cmass = luM^{s)U with U = (p' , iu' , (j)') (11) 

^Their threshold masses are spread out between ~ 1.3 GeV and ~ 2 GeV. The situation is similar for the loops 
generated by the Yang-Mills term. 
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(the ideal fields being supposed real) ancj^ : 



- /iei 



- /ie2 



(with /i = zy^2 ) (12) 



v 



/iei 



where we have defined : 



= (72^^(n+(s)-no(s)) 



< 



^?p/.i.(n+(s) + no(s)) 



(13) 



n. 




In the region where we work - invariant masses bounded essentially by the two-pion thresh- 
old and the </) mass- , the functions e\ and e2 are small and can be treated as perturbation pa- 
rameter^ ; moreover, they are real for real s up to the two-kaon threshold region. In contrast, 
n7r7r(s) is complex starting from the two-pion production threshold and is not expected to be 
small enough to be consistently treated as a perturbation parameter. 

The physical vector meson mass eigenstates are the (s-dependent) eigenvectors of M^(s) 
and their masses are the corresponding eigenvalues, which are also s dependent ! Expressed 
this way, the notion of vector meson mass looks a little bit paradoxical, however, it is not really 
new : writing, as usual, the inverse p dressed propagator Dy{s) = s — rr? — Tlpp{s) can be 
interpreted as stating that the p mass squared is + npp(s) and includes an imaginary part. 
From a physics point of view, what is important is that the pole position associated with the 
p is always a zero of s — — npp(s) located on the unphysical sheet, close to the physical 
regiorlfl. The p pole position has been fitted long ago by [|56ll in e"*"e^ data, and more recent fit 
values can be found in ifSTl ; this piece of information is actually highly model independent, in 
contrast with any other definition [i53il . We shall revisit this issue with our fit results. 

One may wonder about the hermiticity properties of the Lagrangian modified as proposed. 
As below the two-pion threshold, the loops defined above are all real for real s, we still in- 
deed have £(s) = 0{^s), however, above this point, the hermiticity should be redefined as 
£(s) = £l^(s*). This property known as hermitian analyticity [f58l is fulfilled by our modified 
Lagrangian as it is already fulfilled by the loop functions. 

Now, in order to define the physical p, tu, in terms of their ideal partners, one has to 
find the eigenstates of the squared mass matrix given by Eq. (fT2l) . Let us take advantage of 
the smallness of ci and to solve the problem perturbatively in order to avoid dealing with 
untractable expressions. Let us split up the squared mass matrix into two pieces and write it 

''For ease of reading, t\ t2 are not written with their expHcit s dependence which is (or may be) understood 
throughout this paper 

Actually, from their very expressions in terms of kaon and K*K loops, one may expect t\ to be sensitively 
smaller than 62 in absolute magnitude. 

^The upper lip of the physical region- the s > Am^ semi-axis - located on the physical sheet is topologically 
close to the lower lip in the unphysical sheet of the Riemann surface ; in contrast, the lower lip in the physical 
sheet is topologically far from the upper lip in the same sheet. 
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M2 



Mi 



eB with : 



V 













+ €2 







eB 



£2 ei 



ei 



-/iei -/ie2 / 



(14) 

In this splitting up, we have found appropriate to leave a part of the actual perturbation 
inside Mq . In this way, we avoid to some extent the problem of having the unperturbed eigen- 
values degenerated twice or even three times (when zy = 1) for some values of s and some 
zy. However, while assuming that IIt,^{s), ei(s) and £2(5) vanish at origin, one cannot avoid 
to have a twofold degeneracy at s = ; this degeneracy is resolved as soon as s departs from 
zero by an arbitrary small quantity. This issue, which affects strictly the point s = (where the 
exact solution is trivial !), does not raise any problem with our data which are all at s 7^ 0, even 
if close to zero, as the NA7 spacelike form factor data [[59ll . Another solution to this problem 
would be that the HK masses for p/ and cui would be slightly different ; such a mechanism 
remains to be founcj^. 

The unperturbed solution is then trivial, as the eigenvectors are the canonical ideal combi- 
nations of the neutral vector meson fields, with eigenvalues as can be read off the diagonal of 
Mq . Then, one has to solve the following system for the perturbations 5vi and 5Xi : 



Vi = XiVi ,Vi-Vi 



1,2,3) 



(15) 



M\vi + 6vi) = (Ai + 6Xi){vi + 6vi) , with 
for each i = (p, c<j, 0) = (1, 2, 3). The solution can be written : 



Ris) 



Ris) 



Vi ■ 6vi = 



UJ 



(16) 



where (recall Cj = ej(s) are analytic functions of s) : 



( 



R 



1 



pei 



n 



n 



■TTTT \ 



62 



— 62 

1 



;i - zy)m^ + n^^(s) - p2e2 

pea 

1 - zy)m^ + (1 - P^)e2 



pe2 



Pfi 

\ (1 - zy)w? + n^^(s) - y?e2 (1 - zy)w? + (1 - p2)e2 



(17) 



^ A way to get it would have been to use as breaking matrix Xy — Diag(l + £„, 1 + e^, V^y) instead 
of Xy — Diag(l, 1, yiv) while computing £v (see Appendix ICTi. This, actually, generates a mass difference 
between and p^, but the HK mass for the ld meson remains equal to that of the meson. Additionally, the 
coupling constants of the charged and neutral p mesons to a pion pair differ only by terms of order e^y^ and thus 
can be kept equal. 
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The matrix R is orthogonal up to (neglected) second order terms (see Section 6 in [16J) and 
its elements are, actually, meromorphic functions of s ; this, for instance, means that one has 
to check that they do not develop singularities in the region of physical interest for our model. 

On the other hand, one may wonder getting R(s) with no complex conjugation in the 
field transformation (n^^(s) is complex and fulfills the real analyticity condition 11^^ (s) = 
n*^(s*)). This is due to the unitarity condition which writes \Qj : 

R{s + ie) ■ R^{{s + le)*) = R{s + is) ■ R^s - ie) = 1 (18) 

for real s above threshold and e > 0. The real analyticity property fulfilled by the matrix 
function R then gives R*{s — ie) = R{s + ie) and then Eq. (fTS]) becomes : 

R{s + ie)R{s + ie) = 1 (19) 

as can be checked directly with the R matrix above. 

At first order, the corrections for eigenvalues are not changed with respect to their unper- 
turbed values for to and (p , while for p°, the first order correction is such that the eigenvalue is 
restored to + n7r7r(s) + £2 and is formally identical to the mass squarecj^. Therefore, in 
order to deal with the physical eigenstates p°, to and 0, one has to introduce in the Lagrangian 
^ above, the physical fields as defined by Eqs. (fT6l) using Eqs. (fTTI) . For coupling constants, 
one has to perform exactly as explained in Section 6.3 of [[T6l . using the R matrix above and, 
where appropriate, the ideal coupling constants given in Appendix |El 

In order that this solution should be considered valid, one has to check that the non-diagonal 
elements of the matrix R are small compared to 1 in the whole range of application of our 
model. As they depend on fit parameters, this check can only be performed with the fit solution. 

We will not go into more details with expressing the full Lagrangian ([3]) in the basis of 
physical neutral vector meson fields, as formulae become really complicated (even if they can 
be readily and tediously written down). Let us only give the most interesting piece in terms of 
physical fields for illustrative purposes : 



^ag I ^ + ^ag 

P ■ TX a TT = — 

2 2 



, /^ei 

p - — — 1—. + 



(20) 

This clearly shows how kaon loops generate couplings of the physical oj and cf) fields to 
7r^7r+ which vanish (with ci) in the Isospin Symmetry limit. As Y[^^[s) has a large imaginary 
part, it is clear that the phase of the uj coupling compared with p° will be very large at the p 
peak. It should also be mentioned that the couplings shown here (and the matrix elements of R) 
have all a finite limit at s = even if the loops individually vanish at s = as the pseudoscalar 
pairs couple to conserved currents [[SUfTTI . 

The effects of the neglected loops could be briefly mentioned here. The most important 
effect in the expression for R{s) (Eq. (flTl) ) is on the denominators of Ri2{s) and R2i{s) where 
the difference of the anomalous contributions to self-energies for the p and uj mesons will add 
to the present n^^(s) — €2(5) ; this could change a little bit the behaviour near s = where all 
loops tend to zero. 



^For , the mass squared value contains what was named I\'^^{s) and €2 — > 25^^^n_|_/o(s). 
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As stated above, at first order in perturbations, the squared mass eigenvalues are the entries 
in the diagonal of M^(s) given in Eq. (fT2l) . For further use, let us also give the second order 
corrections to the eigenvalues (and thus to the running squared meson masses) : 



el 



2^2 



+ 



/i e 



+ 



n^7r(-s)-e2 (1 - zy)"^^ + (1 - /^^)e2 



(21) 



/i^e2 



[l - zv)m? + (1 - /(i2)e2 



(1 - zv)iv? + n^^(s) 

In the mass range where we are working (from the two-pion threshold to the mass), the 
mass eigenvalues at first order are reai§ for the uj and 0, excluding a width. At second order, 
one observes that the pion loop generates an imaginary part to these mass eigenvalues. Let us 
remind the reader that, as the pole positions are the solutions of s — Aj(s) = 0, one expects to 
find the p pole position close to the value found by [56] . However there is little chance that the 
UJ width happens to carry the correct width value as this should be generated by the anomalous 
pTx loop with additional insertion of a pion loop on the p leg (or simply considering the dressed 
p propagator) or directly through a possible ^ Svr ^ (double) loop. Finally, in the model 
we use, it is only at second order that a difference between neutral and charged p propagators 
(and thus masses) occurs and this is a net (and small) effect of the neutral vector meson mixing. 
This comes in addition to other sources of p° — p"^ mass difference (see footnote (6]). 



6 Renormalization Conditions On Loops 

With the approximations we choosed (neglecting the anomalous loop contributions), the 
loop expressions needed in order to construct the pion form factor F^(s) are only the vrvr and 
KK loops. They can be computed by means of Dispersion Relations ['52', 9] and can be derived 
without explicit integration, relying only on properties of analytic functions, especially the 
uniqueness property of the analytic continuation (see|^ Appendix A in 

From general principles, any loop n(s) is a so-called real analytic function (see the Section 
just above), the imaginary part of which is calculable using the Cutkosky rules, or in the simple 
case of single loops, using the partial width of the process V PP' . Indeed : 

Imn(s) = -v^r(s) , s > So (22) 

where s is the pair invariant mass squared, sq the threshold mass squared of the pair and V 
the partial width of the decay. With this at hand, the full loop is the solution of the integral 
equation : 

s" /""^ Imn(2; 



n(s) = P„_i(s) + - / — "-±—,dz (23) 

TT Jso z^(z — s + le) 



^Actually, near the accepted (j> mass our perturbation parameters start to carry a tiny imaginary part. 

^In this reference, the loop expressions for equal mass pseudoscalar meson pairs and vector-peudoscalar pairs 
are already given and are correct ; the function given for unequal mass pseudoscalar meson pairs is not correct 
as the contribution of the gauge term has been omitted ; we apologize for this error and correct for in the present 
paper (see Appendix|F|. 
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where P„_i(s) is a polynomial of degree n — 1 with real coefficients and the integral runs 
over the right-hand cut (the physical region). These coefficients should be fixed by means of 
condition to be fulfilled by n(s), the so-called subtraction polynomial, which are nothing but 
renormalization conditions chosen externally and depending on the problem at hand. 

A priori, the number of subtractions, i.e. the number of conditions can be arbitrary, however, 
in order that the integral in Eq. (|23l) converges, there is a minimal number of subtractions to 
perform : For PP' and VV loops n > 2, for VP loops n > 3. 

In the most general form of the pion form factor following from the Lagrangian (|3]), and 
using the modified one-loop mass term (flOl) the relevant basic loops are only the pion loops and 
the kaon loops. These imply that at least n = 2 ; however, the very existence of VP neglected 
loops implies that we are still minimally subtracting with using n = 3 for all loop functions 
in this paper. As discussed in Section |4l in this way, the subtraction polynomials carry some 
(unknown) information on the anomalous loop contribution. 

Additionnally, we request all polynomials P„_i(s) to fulfill P„„i(0) = reflecting this way 
current conservation Il52l [TTl [6T1 [STTl when needed and an appropriate constraint otherwise^. 

In usual approaches [t6n [571 l52l . the renormalization conditions are defined from start and, 
then, one leaves free other parameters like meson mass and width in order to accomodate the 
experimental data. As already done in lISTIl . we proceed in the opposite way : as masses and 
couplings are fixed consistently in our Lagrangian, we leave free the subtraction polynomials 
in the loops n^7r(s) , ei(s) and £2(5). This way allows the full data set to contribute to fixing the 
subtraction constants. The basic (pion and kaon) loop expressions are given in Appendix E and 
are used only subtracted once (in order that they vanish at s = 0) ; they are supplemented with 
second degree polynomials vanishing at s = and having coefficients to be fixed by fitting the 
data. 



7 The Model Pion Form Factor In r Decay 

Introducing pion and kaon loop effects gives the p"^ a self-mass, but, nevertheless, the p"^ 
fields remain mass eigenstates. To stay consistent with using p self-mass, one has also to 
account for loops in the W — p transition amplitude. In r decay the relevant loop effects, while 
neglecting anomalous VP loops, are the vr 
modifies Eq. (01) to : 



±^0 



and K^K^ contributions. Accounting for this 



a, 
2' 



F^p 9pinT 



1 



(24) 



with : 



m 



Kp{s) 



(25) 



where = ag'^ fl and the charged p self-mass n^^(s) has been defined in Section |4]by Eq. 
^ and used in Eq. (flOl) (recall we neglect VP loops). One should note that Eqs. (|24l) and (l25l) 



'°In tliis case, tlie constant term in the p propagators is the squared (HK) mass occuring in the Lagrangian with 
no modification. 
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are not affected by any breaking mechanism. The diagrams contributing to the pion form factor 
in T decays are sketched in Figure [U 

Let us denote for a moment the pion and kaon amputated (i.e. computed with unit coupling 
constants) loops by £,r(s) and ixis), assuming they are already subtracted once in order that 
they identically vanish at s = (see AppendixlF]). The W — p transition amplitude and the 
self-energy occuring in the pion form factor have the following expressions in terms of pion 
and kaon amputated loops : 



^W{S) — QpTTTT 



+ Pwis 



24 



(26) 



+ Pais) 



where (yf^TTTr = afl'/2 and Pwis) and Pp{s) being subtraction polynomials with real coefficients 
to be fixed by external renormalization. As emphasized in [51] and [61 J . the polynomials 
Pw{s) and Pp{s) can be chosen independent. Indeed, Im Ilw{s) and Im Il'pp{s) are even not 
proportional as soon as SU(3) is broken (za 7^ 1) ; moreover, the transition amplitude nvi/(s) 
is non zero even if a = 2 as soon as SU(3) symmetry is broken. We choose to constrain Pwis) 
and Pp{s) to be second degree and vanishing at s = 0, as discussed in Section[6l 

For the sake of simplicity, we have also chosen to use the pseudoscalar meson loops assum- 
ing m^± = m^o and m^± = uiko after having checked that this is numerically armless while 
dealing with all form factor data. Under this approximatior0, we have Ii'pp{s) = Ilpp{s) (see 
Eq. ([8])) and all symmetry breaking effects due to the pion mass difference are concentrated in 
the phase space factors (see Eqs. dSHV])) for cross sections and partial widths where these have 
a sizable effect. 

Therefore, one can rewrite Eqs. (|26l ) under a form more appropriate for our fitting proce- 
dure : 



pp\ 



9pTT 

n 



'1 



-)n^i 



s) + 



(27) 



+ e2 s 



where €2(5) 
defined : 



already defined in Eqs. ([T3l ). carries its own subtraction polynomial, and having 

) + Qwis) 



7 pTTTT 



[■V = 9. 



pTTTT 



{s) + Qp{s) 



(28) 



where Qwi-s) and Qp{s) are second degree polynomial^ (with real coefficients to be fitted) 
and vanishing at origin. Possible correlations among them, if any, would be an outcome of the 
fit procedure and can be detected from inspecting the fit error covariance matrix. Finally, one 
can check that the condition F^{0) = 1 is automatically fulfilled 

Before turning to the pion form factor in e+e annihilations, let us also remind the reader 
that the r partial width expression in Eqs. ([5]) has to be further corrected for isospin breaking 



' ' This implies that the and widths do not significantly differ. This statement is supported by the various 
experimental data collected in f5T| . 

'^We recall here, that these polynomials may account for the neglected anomalous loop effects not introduced 
explicitly. 
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effects by multiplying it b)o Sew = 1 0232 which accounts for short range radiative correc- 
tions ll37l . Long range radiative corrections have been derived in JSHl [39l |40l| and come as a 
further factor Gem{s) ; another estimate taking into account additional Feynman diagrams can 
be found ll4Tll42ll and a corresponding numerical parametrization of Gem{s) has been provided 
in ||43l . This means that in all our fits we perform the substitution : 



FKs) 



SewGem{s)FI{s) 



(29) 



which, therefore, accounts for all reported Isospin Symmetry breaking effects specific to the 
r sector. Another isospin breaking effect might have to be considered, namely a p° — p"*^ 
mass difference. This can be generated, for instance, by means of the mechanism sketched in 
footnote [6l It can be shown that this turns out to modify Eqs. (|25l) to : 



r f; = r - Yiw{s) 



DM 



s — m^ — dm^ — TI'pM 



(30) 



5m? 



where Sni^ is left free. The modified Eqs. xd® allows Eq. ^ to still fulfill F;(0) = 1 
automaticall>0. 



8 The Model Pion Form Factor In e+e Annihilations 

In r decay, the pion form factor, as just seen, is free from any vector meson mixing effect. 
Instead, the pion form factor F^{s) expression is sharply influenced by the vector meson mix- 
ing mechanism constructed explicitly in Section [5] which leads us to make the transformation 
from ideal to physical vector meson fields. After this tranformation, we get from our effective 
Lagrangian the diagrams shown in Figure [Hand the corresponding expression : 



(31) 



where Dp(s) (seeEq. (l25l) for its charged partner). D^(s) and-D<^(s) are the inverse propagators 
of the corresponding (physical) vector mesons. We have now : 

D,{s) =s-m'- Iipp{s) (32) 

(recall that our assumptions on pseudoscalar meson masses implies Tipp{s) = H'^p^s), which 
reduces the number of free parameters in our model). The vector meson couplings to a pion 



'^Actually, this numerical value has been derived for the pion final state ; in practical applications, it is usually 
assumed that this value holds also for the p final state - see for instance p4 , 35, 361. 

''^ Anticipating somewhat the fit results, a possible Sm'^ can be detected on ALEPH data ITl | (not on CLEO data 
||33l ) and amounts to ~ -0.25 GeV^. /^(~ 0.7 GeV^) is then increased by 0.5 lO'^ GeV^ a quite negligible 
quantity. 
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pair after symmetry breaking are : 

9p7V7T 
9 LOU TV 



2 



ag 



2 



ag 



£2 

/iei 



(33) 



2(1- zy)m2 + n 



where n^^(s) has been defined in the previous Section. One should note that the quantity 
named Yi-„T,{s) in the definition of the matrix R{s) (see Section[5]) coincides with the presently 
defined n^^(s). 

The quantities Fy can be written : 



n 



V71 



(34) 



Collecting the various couplings of the ideal fields suitably weighted by elements of the 
matrix transformation R{s) (see Eq. (flTl)). we get : 



/; = agfl 



1 

1 + 77 



ei 



3 n^7r(s) - 62 

1 H 

3 n^7r(s) — 62 

3 



1 

3(T 



TTTT (5) 



1 



/i^e2 



3 (1 - 2y)m2 + (1 - /i2)e2 



(35) 



/i6i 



/i 



£2 



(1 - 2y)m2 + n 



TTTT (s) 



3 (1 - 2y)m2 + (1 



/^^)e2. 



and, keeping only the leading (first order) terms, the loop corrections Uyyis) (see the defini- 
tions in Eqs. (fT3l) . and the expression for /i in Eqs. (fT2l) ) are : 



npo^(s) 



n^^(s) 



n^^(s) 



2 gp-KTT 

-(!-- 



+ [ZA 



a ,.61 + 62 
7T -^)- 



£2 - ei 



n;[.(3) 



2'n^^(s)-62 gpn 



+ {zA---b) 



QpTTTV 

ei + 62 



(S) - /i^62 gp^ 



, £2 - 61 

— 

gpniT 

f - 

2 gp^ 



£2 



£2 



(36) 

The first term for each transition loop is the pion loop contribution while the others are 
resp. the charged and neutral kaon loops. Of course, the functions occuring there are the 
same as for F^. We have denoted by n^^(s) the transition amplitude for 7 — , which is 
in correspondence with the W — transition amplitude introduced in the previous Section 
(see Eq. (|28l)). A priori, the subtraction polynomials of n^^(s) and n)j^(s) might be slightly 
different. However, in an attempt to reduce further the number of free parameters of the model, 
we assume that they coincide, which turns out to identify the amputated W — and 7 — 
transition amplitude^ We shall see that this assumption is well accepted by the data and, 
moreover, make clearer the switching to the r form factor expression. 



'^This is a nothing but a strong CVC assumption. 
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In addition to the explicit dependence of our model on the HLS basic parameters a, g, and 
on the breaking parameters x, za, zv,zt and Sm?, there is a further dependence on subtraction 
parameters hidden inside 11^^ (s), IlYJ'^{s), ei(s) and £2(5). Isospin symmetry breaking is 
reflected in having a non-zero ei(s) function. We have : 

' ^Tis) =Qo{s)+iAs) 

(37) 

e,{s) =P_{s)+iK±{s)-eMs) 

62(5) = P+is) + iK±is) + iKoi-s) 

where i^is), iK±{s) and iK"{s) are now the non-amputated ix^ix^ , K^K^ and K^K^ loops, 
subtracted in order that these loops vanish at the origin. The parameter polynomials Q^i^s), 
Pq{s), P-{s) and P+{s) are chosen to be second degree with zero constant terms in order to 
stay consistent with the Node theorem ifTTllSlll . 

One can check that -F^(O) = 1 + 0{el), which could have been expected from having 
neglected terms of order greater than 1 in our diagonalization procedure^. 

As the form factor data collected at the (f) are not currently available, the last term in Eq. (|3T| ) 
could have been removed. However, in order to account for tails effects, we preferred keeping 
it and use a fixed width Breit-Wigner expression incorporating the Particle Data Group mass 
and width recommended values [54]. Due to the narrowness of the u mass distribution, we also 
have replaced in our fits the to propagator by a fixed width Breit-Wigner constructed using the 
recommended mass and width from [54J. 

Let us also recall that Eq. (|3TI) is our formula for the pion form factor F^{s) in both the 
spacelike and timelike regions. Indeed, for consistency, we will not remove the to and meson 
contributions while going to negative s. 

Using the first order correction to the p mass eigenvalue, the inverse propagator could 
have been written -Dpo(s) = s — Xp{s), as the leading order squared mass eigenvalue is : 

A,(s)=m2 + ne,(s) + e2(s) (38) 
As far as e+e" data are concerned, we shall modify the eigenvalue expression to : 

X,is) = m' + UP^is) + e^is) + S^X^is) (39) 

by adding the second order correction given in Eqs. (|2TI) . This does not add any more freedom 
in the model, but rather allows some check of the diagonalization method. 

Therefore, the difference between F^{s) and F^(s) is solely concentrated in the coupling 
changes from ideal to physical fields given by the varying matrix R(s) (see Eq.(fT71)) which 
only affects -F^(s). Stated otherwise, modifying the function F^{s) in order to incorporate 
isospin breaking effects is strictly equivalent to using our expression for F^(s) directly, the 
factor SewGem{s) being removed and Srn^ being made identically zero. 

Actually, it depends on the (1,2) rotation matrix element : 1 + 0{[Ri2{s — 0)]^) and, numerically, the 
neglected term is ~ 1.5 10^^. 
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9 The Photon Vacuum Polarization 



The raw data on the pion form factor F^{s) should be "undressed" by unfolding the con- 
tributions due to radiative corrections and to the photon vacuum polarization (VP) before any 
comparison with r data (we refer the reader to [62l |63l for a comprehensive analysis of these 
factors and for previous references). Quite generally, available experimental data on -F^(s) 
have already been unfolded from radiative corrections [|27l |28l |29l |30l EH [63|. All the data 
sets just referred to are not unfolded from photon vacuum polarization (VP) effects, except for 
KLOE data 1(641 . Therefore, one has to account for VP effects by including the corresponding 
factor when comparing a pion form factor model with experimental data. Traditionally (see for 
instance [65] and references quoted therein), this results in the change^ : 

F:{s) ^ (1 - Uypis))F:is) (40) 

when comparing with most data sets. 

The VP function Ilvp{s) contains two parts. The first one is the sum of the leptonic loops 
which can be computed in closed form at leading order (In Appendix |Fl we recall the explicit 
form at order a and give its expression along the real s axis). The second part is the one 
particle irreducible hadronic contribution to the photon self-energy which is derived by means 
of a dispersion relation ; at low energy, where non-perturbative effects are dominant, this is 
estimated using the experimentally mesured e+e~ cross section (see, for instance, If65l [6611671 
l62l ), while the high energy tail is calculated using perturbative QCD. 

For our purpose, we use the sum of the leptonic VP as given in Appendix |F| for e+e", /i+yU^ 
and r+r together with a numerical parametrization of the hadronic contribution. From the 
2-pion threshold to the (f) mass, we benefited from a parametrization0 provided by M. Davier 
[|68l . Below the 2-pion threshold and down to s = —0.25 GeV^, we use instead a (real valued) 
parametrization provided by H. Burkhardt [|69l . 

10 Decay Widths Of Light Mesons 

In order to compute decay widths and fit data, one has to define the couplings allowing the 
decays of the light mesons involved. For the two-photon decays of the i] and r]' mesons, as well 
as for the radiative decays of the and A'*'s mesons, the couplings defined after SU(3) and 
Nonet Symmetry breaking (see Eqs. (11031) and (11061) ) are the couplings coming directly in the 
decay widths formulae (see Subsection IE. 31 ) and do not depend on further Isospin Symmetry 
breaking effects than mass values in phase space factors. 

The isospin breaking procedure we presented plays only for the p°, cu and mesons. In 
order to compute the leptonic decays of these, one has to use the full couplings Fy^s) as 
given by Eqs. (|35l) and (l36l) in Eq. (II 101) and computed at the appropriate vector meson 
masses Fy(my). As the loop functions are slowly varying, one can choose niy as the Higgs- 
Kibble masses occuring in the Lagrangian (see Eq.([T])), which moreover simplifies the fitting 
procedure. 

If one denotes by I](s) the photon self-mass, the inverse photon propagatoris given by _D;:^^(s) = s — S(s) = 
s(l — I](s)/s). Therefore, compared with traditional notation, we have Hyp (s) = — I](s)/s. This is not a problem 
but should be kept in mind. 

'^This parametrization neglects the (very small) imaginary part of the hadronic VP contribution. 
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For the other (radiative decay) coupling constants one has to combine the ideal coupling 
constants (given in AppendixlE]) using the transformation R(s) to derive the physical coupling 
constants, as was described in [fT6l ; the context, compared to [fT6ll . slightly differs due to the 
fact that, now, the mixing parameters are functions to be computed at the appropriate my values 
in order to go to the vector meson mass shell. 

Traditionally, the p° is decoupled from mixing and treated as the and mixing effects are 
only considered in the (cu — 0) sector. Additionally, it is usual to treat the (u — (p) mixing angle 
as a constant to be fit (see [H |5l HI [H [161 and the references quoted therein). The approach in 
the present study is different : One considers a full (p° — — 0) mixing scheme (as in [[161 ). 
however - for the first time - the mixing parameters are functionally related and the same 
functions have to be computed at each vector meson mass. For instance, the (cu — 0) mixing 
"angle" has not the same numerical value at the lu mass and at the mass. This only reflects 
that the mixing is actually invariant mass dependent. When, as for — tttt [|50l and — tttt 
[[70l [TTl data exist on the phase of the coupling constant, these phases can be introduced in the 
fit with the same functions, the modulii of which determine the branching fractions. 

1 1 The Full Set of Data Submitted To Fit 

In order to work out the model presented in the Sections above, we use several kinds of 
data sets. In this Section, we list them and give some details on the way they are dealt with in 
our fit procedure. 

11.1 Partial Width Decays 

As a general statement, the decay data submitted to fits have been chosen as the so-called 
"fit" values recommended by the Particle Data Group (PDG) in the latest (2006) issue of the 
Review of Particle Properties [|54l . 

This covers, with no exception, the leptonic decay widths of the p°, cu and mesons, the 
two-photon decay widths of the r] and 7]' mesons and the tt+tt^ partial width of the meson. 
There are two measurements of the phase of the coupling constant g^.,r+Tv- reported in the 
literature ; the older one ^is^^ = -42° ± 13° and more recently nB^^ = -34° ± 4° ± 3°. 
Summing up in quadrature the errors, we choose as reference value in our fits tptf, = —34° ± 5°. 
We could have chosen to include in our fits the 2-photon decay width of the 7r° ; however, we 
preferred replacing this piece of information by the pion decay constant value f^, = 92.42 MeV 
and did not let it vary, as this is supposed to carry a very small erroi0 [|54l . 

The RPP pieces of information [|54l concerning the p and cu decay width to tt+tt^ and the 
partial width p° — > e+e^ are not considered as data to be submitted to fits, as they have all been 
extracted from fitting the same pion form factor timelike data which we are included in our 
fit procedure (see the Subsection below) ; this does not prevent us from comparing our results 
to the RPP available information. This is also true for the relative phase of the couplings 
9ujTT+TT- to gpTT+TT- (^c so-callcd Orsay phase) which has been measured [|50l with the result 

= 104.7° ±4.1°. 

'^Possible Isospin Symmetry breaking effects might have to be considered. 
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Instead, it is quite legitimate to include the uj e+e^ partial width in our fit procedure 
as, even if this mode could have been marginally influenced by the pion form factor data, it is 
mostly extracted from e+e^ tt+tt^tt'^ data [[54l . As the pion form factor spectrum around 
the mass is not currently available, the — > e+e~ partial width is quite legitimately included 
in our fit data set. 

In the present work, as in our previous works on the same subject [[8l|9l|T6l, we do not intend 
to use the decay widths K* Kn. Actually, as for the width p — »• tttt which is inherently 
fitted with the pion form factor, the choice of the mass value for a very broad object makes the 
extraction of coupling constants a delicate matter. It should be more appropriately discussed 
with the Ktt form factor in r decays when the corresponding data will become available. 

The data on the two kaon partial widths of the meson are also left outside our fit proce- 
dure. In a previous work of some of us [fT6l . as in other works [|72l the issue of accomodating 
the K^K^ partial width was raised. A recent work [|73l claims that the ratio of these par- 
tial widths can be accomodated by introducing corrections to the decay widths which increase 
both partial widths as derived from the matrix elements of the transitions. As then, the problem 
may affect both the — > K^K^ and the — > K^K^ partial widths, we have preferred leaving 
both modes outside the fit procedure. We will discuss this issue below in a devoted Subsection. 

We also use all radiative decay partial widths of light flavor mesons of the form V P7 
or P ^ V'j. As a general rule, we chose as reference data the "fit" values recommended by 
the PDG as given in the latest RPP issue [54]. There are two exceptions to this statement : the 
partial widths for ^ 777 and to tt^'J. 

Indeed, as already noted in [fT6ll . there is some difficulty in accomodating the present u 
r]'j "fit" branching fraction ((4.9 ± 0.5) 10~^) while the so-called "average" value [54] ((6.3 ± 
1.3) 10~^) is much better accepted by our model fitting. 

On the other hand, the new "fit" value for the branching fraction cu — > 7r°7 ((8.90 ± 
0.25) 10^^) is also hard to accomodate in our model. We show that the previous PDG "fit" 
value ([8.50 ± 0.50] 10^^) seems in better consistency with the rest of the data we handle ; we 
also checked that the value ([8.39±0.25] 10^^) produced by a fit [|74l performed in a completely 
different context was also well accepted, pointing to a possible overestimate of the central value 
for this mode0. The questions raised by the values of these two decay widths will be discussed 
at the appropriate place below. 

Finally, we also introduce in the fit procedure the ratio of the kaon to pion decay constants 
as they are reported in the latest RPP [|54l . This actually coincides with our SU(3) breaking 
parameter (z^ = [/k//^]^). 



11.2 Timelike Form Factor Data in e+e Annihilations 

Beside the decay data listed just above, we have included in our fit all data on the pion form 
factor collected in e+e^ annihilations by the OLYA and CMD Collaborations as tabulated in 
iTTSl and the DM1 data [76] collected at AGO (Orsay). These data will be referred to globally as 
"old timelike data". When included into a expression, systematic errors have to be combined 

^"Looking at |54|, the role of the data and analyses for the e+e^ — > tt'^j process itself to get the "fit" value for 
the partial width w ^ 7r'^7 is unclear. 
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with the published statistical errors ; they have been first added in quadrature to the statistical 
errors for OLYA data (4%) and CMD (2%) following expert advice [|77|| . However, for sake of 
consistency with the new data discussed just below, we preferred extracting the correlated part 
of the systematic errors, estimated fTTl to 1% and have performed the same treatment as for 
the new data (see just below). The accuracy of the DM1 data making the influence of this data 
set marginal, we did not add any further contribution to the published errors. We only use the 
data points located below the cj) meson mass in order to avoid being sensitive to higher mass 
vector mesons, not included in the present model. 

Four additional data sets have been collected later at Novosibirsk on the VEPP2M ring. The 
first one, covering the region from about 600 to 960 MeV, collected by the CMD2 collaboration 
|[78l and recently corrected [27J, is claimed to have the lowest systematic error (0.6%) ever 
reached in this field. 

CMD2 has collected in 1998 and recently published two additional data sets, one ll28l 
covering the energy region from 600 to 970 MeV is claimed to reach a systematic error of 
0.8%, and a second set [29] covering the threshold region (from 370 to 520 MeV) with an 
estimated systematic error of 0.7%. On the other hand, the SND collaboration has recently 
published ll30l a new data set covering the invariant mass region from 370 to 970 MeV with a 
systematic error of 1.3% over the whole data set except for the very low mass region where the 
(first) 2 points carry a systematic error of 3.2%. 

Concerning these four data sets (which will be referred to globally as "new timelike data"), 
we could, as per usual, add in quadrature the systematic and statistical errors and then get a 
diagonal error matrix which can be used in fits in a trivial way. 

However, an important part of the systematic uncertainties for these data sets is expected 
to be a common global scale uncertainty IfTTl which has been estimated to 0.4% and generates 
bin to bin correlated errors. In principle, one should take the latter information into account in 
fits ; this implies dealing with systematic and statistical errors in a way slightly more elaborate 
than simply adding in quadrature statistical and systematic errors. 

Firstly, the (bin per bin) uncorrelated part of the systematic error is derived by subtracting 
in quadrature 0.4% from the already quoted systematic errors. This uncorrelated part of the 
systematic error {i.e. ^Jcrgyst. ~ (0.4%)^, depending on the data set considered) can certainly 
be added in quadrature to the statistical error bin per bin, giving a combined standard deviation 
named for the measurement in the energy bin i ; the ai are uncorrelated errors and define 
a diagonal error matrix. The question then becomes how to redefine the full covariance matrix 
for each experiment, being understood that the quantity to be compared with the theoretical 
pion form factor for each energy bin i is related with the measurement rrii by : 

/f ^ = (1 + 5A)m, (41) 

where 5\ is considered a gaussian random variable with zero mean and standard deviation 
A = 0.4 X 10^^. With this assumption it is possible to model reasonably well the covariance 
matrix, which is no longer diagonal. 

Secondly, one has to treat these correlations. The quoted correlated systematic error is a 
conservative estimate of the accuracy of radiative corrections performed on the four data sets 
using the same Monte Carlo generator flT]. Therefore, the fit parameter introduced in order 
to optimize the absolute scale should be the same for all data sets ; in statistical terms, this fit 
parameter value can be considered as only one sampling of the gaussian random variable 5\ 
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defined just above and should be valid for all data collected by the CMD2 and SND Collabo- 
rations. 

If the correlated part of the systematic error was strictly zero, the error covariance matrix 
for each data set would simply be given by : 

V^, = (42) 

where i and j label energy bins in the data set. In the case of existing correlations, having 
defined Oi as the sum in quadrature of the statistical error and the uncorrected systematic error 
in the z*'^ bin , the error covariance matrix elements can be written : 

Vi,=Y.MiuWuMi^ , (43) 

k,l 

where : 

Mij = ai6ij and , Wij = 6ij + X^CiCj (44) 

A = 0.4 10"^ being the standard deviation of the correlated error function and the vector e 
being defined by its components on the various energy bins i as the ratio of the corresponding 
measurement to its uncorrelated error : 

Ci = rrii/ai , Vz G [1, ■ ■■n.measur] (45) 

However, what is relevant for fitting, is not so much the covariance matrix Eq. (I44l) as 
its inverse. It happens that the matrix W can be inverted in closed form : 

W^^= S,j-fi^eie, and /x' = ^ _^ j^^—r. ^2 (^6) 
and the full error covariance matrix is also inverted in closed form : 

V^' = J2Mr,'W,fM[:' . (47) 

k,l 

This is, together with the measured values, the main ingredient of the calculation which 
will be performed with the four new timelike data sets. Finally, while fitting the new data, a 
term has to be added to the \ naming A fu the fit parameter for the global scale common to all 
the new Novosibirsk data sets and Ae^p = 0.4 x 10~^ the scale uncertainty on the measured form 
factor estimated by the experiments, this additional contribution to the is simply [A fu / Ae^p]^- 

Mutatis mutandis, the same method has been applied to the old Novosibirsk data sets using 
another global scaling factor Aj^^ with Ag^.p = 1.0 10^^, as recommended by informed people 
IfTTl . and the same procedure to construct the final inverse covariance matrix to be used in fits. 

A new data set has been recently collected by the KLOE collaboration [|64ll using the Ra- 
diative Return Method. Existing analyses (see, for instance, the short account in ||26l ) however 
report a disagreement between KLOE data and the recently collected data sets at Novosibirsk 
due to some systematic effect presently not understood. A recent study of a parametrization of 
the pion factor [79] argues about a possible systematic energy shift in the data which would be 
detected by fitting the to mass. In view of this unclear situation, we have found it appropriate 
to postpone including the existing KLOE set among our fitting data samples. 
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11.3 Spacelike Pion Form Factor Data 



In order to further constrain the pion form factor in the timelike region, information on 
the close spacelike region is valuable. Reliable data on the pion form factor in the negative s 
region are somewhat old [|59l[80ll . The Fermilab data set [80] consists of 14 measurements of 
|F^(s)|2 between s = -0.039 GeVVc and s = -0.092 GeVVc with 2^7% statistical error ; 
an estimated systematic error (overall normalization) of 1 % is provided in [[80ll . The NA7 data 
cover the region between s = —0.015 GeV^/c and s = —0.253 GeV^ with 45 measurement 
points and an overall statistical precision better than those of the Fermilab data. However, NA7 
data are also claimed to undergo an overall scale error of 0.9 % rms. 

One will use these two data sets and treat these correlated systematic errors exactly as 
explained above for the timelike pion form factor data. 

Data have more recently been collected at the Jefferson Accelerator Facility [|8T| and reana- 
lyzed in order to optimize the extraction of the pion form factor data in the region for s between 
—0.60 and —1.60 GeV^ with a quoted uncertainty of about 10%. No precise information about 
the correlated-uncorrelated sharing of the systematic error is reported. Including these data in- 
volves some more studies and modelling which goes beyond the main task of the present work, 
namely, to check the consistency of e+e~ and r data. 

11.4 Phase Pion Form Factor Data 

There are several data sets available which provide measurements of the isospin 1 part of 
the vrvr amplitude phase shift. The most precise set is the CERN/Munich one [[82l . but the 
older Fermilab data set [,8 3 J is still useful. However, systematic errors here are not completely 
controlled. Moreover, as we neglect vertex corrections at the tttt vertex and t-channel resonance 
exchanges which may carry some unknown imaginary part, one cannot draw firm conclusions 
when comparing the phase information of our pion form factor with phase shift data. Therefore, 
we have left these data outside our fitting procedure and limited ourselves to simply compare 
graphically with the phase of our pion form factor. 

11.5 Pion Form Factor Data From r decays 

There are presently three available data sets concerning the pion form factor. These have 
been collected at LEP by ALEPH fi311 and OPAL [32] Collaborations and at much lower energy 
by the CLEO Collaboration [33l|. 

The data provided by the ALEPH Collaboration [f3T| include the covariance matrices for 
statistical and systematic errors which should be added before inversion in order to be used in 
a minimization. There is some disagreement between ALEPH ||3T| and OPAL data ||32l 
which has led most works to discard this data set ; we shall do likewise. 

The CLEO data [33] on the pion form factor are also provided with their full error matrix, 
but one that accounts for statistical errors only. Statistical errors dominate most of the system- 
atic uncertainties except for those contributing to the absolute energy scale for determining ^/s 
lf84l . These were quantified by CLEO as a systematic uncertainty on the value of the mass 
obtained in their fits to form factor models, estimated to be 0.9 MeV. This error, not accounted 
for by the CLEO error covariance matrix, is a systematic error which correlates the various bin 
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energy values. In contrast with the Novosibirsk data, it is not easy to rigorously account for 
this correlated systematic erro As an approximation we allow the central bin y/s value to 
vary by some e MeV and add to the CLEO data x^- This approach provides a simple 

and reasonable way to deal with the data and errors If84l . 

In order to stay consistent with our dealing with e+e" data, we have limited our fitting 
range to the (f) mass and then removed all points above s = 0.9 GeV^. This leaves us with 
33 measurements from ALEPH and 25 from CLEO, largely unaffected by higher mass vector 
meson effects, as will be checked. 

12 The Main Global Fit To The Data Sets 
12.1 General Comments About The Fits 

Our global model has seven parameters carrying an obvious physical meaning : 

• The universal vector coupling g, 

• the SU(3) breaking parameter za (expected to coincide with [/a'/Ztt]^ within errors), 

• the Nonet Symmetry breaking parameter x, 

• the basic HLS parameter a (expected close to 2), 

• the parameter zy which mostly governs the mass difference between the — to system 
and the meson but also plays a role in some coupling constants, 

• zt which affects only the K* radiative decay sector in the data used, 

• and, finally, the p° — squared mass shift Sni^. 

These have been already fitted in isolation in related previous works (81 |9l [161 US [851 
and we expect to find fit values close to the already published ones. Within our approach, the 
pseudoscalar mixing angle is not free, but is derived from the previous parameters using Eq. 
(11081) and is expected close to —10.5° degrees from previous fits [[8l[T6l. This has been found in 
perfect agreement with the two-angle formulation [11] expressed in the framework of Extended 
Chiral Perturbation Theory [13l[l4l. 

Beside these parameters carrying a clear physical meaning, one has the subtraction polyno- 
mial of the pion loop (mostly associated with the p meson self-energy), assumed to be written 
cis + C2S^ with ci and C2 to be fitted. Two additional subtraction polynomials carrying the same 
form and associated with the difference (ei(s)) and the sum (£2(5)) of the K^K^ and K^l^ 
loops introduce 4 more parameter^ to be fit. Finally, two more subtraction parameters come 
from the specific subtraction of the 7p (or W p^) transition amplitude. We thus end up with 
15 parameter^ for a number of data of 344 (18 decay modes, 127 data points from the new 

^'it affects the position of the measurement, not the measured value itself. 

^^We note that we approximate the K^K'^ loop by the average value of the K^K^ and K'^K loops in order 
to limit the number of free parameters. 

^^It will be emphasized later on that one among these subtraction parameters does not influence the fit and can 
be safely fixed to zero. 



26 



timelike pion form factor data, 82 from the old timelike pion form factor data, 59 data points 
in the spacelike region, 33 data points coming from ALEPH data and 25 from CLEO). 

In addition to these parameters which define our model, we have to account for correlated 
systematic errors in several experiments by fitting the corresponding scale factors and using 
the experimental pieces of information as constraints. These additional degrees of freedom are 
therefore exactly compensated in number by the constraints. This covers the global scale factor 
of the former Novosibirsk experiments as reported in [75] (estimated to 1.0% r.m.s.), the global 
scale factor of the new Novosibirsk experiments as reported in [|27l l28l |29l [30ll (estimated to 
0.4% r.m.s), the scale factor for the NA7 ll59ll and Fermilab [80] data (estimated respectively 
to 0.9 % and 1.0 % r.m.s.). Finally, the CLEO data set is expected to carry a systematic energy 
shift which will be fitted and is expected [84] of the order 0.9 MeV. 

We have performed various kinds of fits. In all of them, as detailed above, we have 
introduced all usual symmetry breaking effects in the value of meson masses, the prominent 
effects oi — uj — (j) mixing (for e+e~ data) and the long- and short-range [|37ll IB correction 
factors (for the r spectra). We have observed that the two proposed ways to account for long 
range corrections by either of ll38l [39l |40I and [411. |42[ l43l approaches provide quite similar 
effects and that, on the basis of probabilities, the difference was never observed significant in 
any of our fits. For definiteness, we choose to use the function of [[38l [39l |40|| for all results 
presented here. 

On the other hand, it was useful to check the effect of excluding the photon vacuum polar- 
isation (VP), by fixing the corresponding factor to 1. We also found it of interest to perform 
fits by excluding either the r data or the spacelike data ; this gives information on the effect of 
these on the global fit quality and on the stability of the fit parameter values. Finally, it has also 
been of interest to check the mass shift effect between the p° and the mesons, by fixing the 
corresponding parameter to zero. The results summarizing the statistical qualities are gathered 
in Table [Hand the fit parameter values can be found in the appropriate data column of Tables |2] 
and [31 

12.2 Discussion Of The Fit Information 

Table [T] reports the statistical information about our fits under various conditions. As a 
general statement, the fit quality is always either reasonable or very good as clear from the last 
row in Table [U 

As a first remark, neglecting to account for photon vacuum polarization effects does not 
end up with a dramatic failure ; however, there is a general improvement while introducing 
the corresponding function. The negligible degradation observed for CLEO data is entirely 
produced by the value found for the CLEO mass shift parameter (which contributes to the x^, 
as explained above) will be discussed below. This statement clearly follows from comparing 
the two data columns named "full fit" in Table [IJ 

The gain in y^, while including the photon VP, is 13.5 units without any additional param- 
eter freedom and, in terms of fit probability, one wins 20 %. Therefore, one may conclude 
that the data description prefers including explicitly the photon VP while fitting the e^e^ data. 
Under realistic conditions, the fit probability is then always of the order 75 % or better. 

The fits have always been performed using the package MINUIT |[86ll and the errors quoted 
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Data Set 


Without VP 


With Vacuum poliarisation (VP) 


H (data + conditions) 


Full Fit 


Full Fit 


No r 


No Spacelike 


Nop mass shift 


Decays (18+1) 


11.46 


11.13 


11.52 


11.48 


11.25 


New 
Timelike (127+1) 


132.81 


128.10 


122.02 


125.76 


132.23 


Old 

Timelike (82+1) 


62.22 


59.05 


54.68 


55.20 


60.15 


SDacelike('59+2') 


68.53 


65.70 


55.20 


89.82/(59) 


65.13 


T ALEPH (33) 


27.06 


23.86 


42.27/(33) 


20.80 


24.48 


r CLEO (25+1) 


25.53 


26.06 


26.16/(25) 


29.72 


28.55 


Probability 


327.40/331 
54.6 % 


313.83/331 
74.3 % 


257.73/274 
75.2% 


238.81/272 
92.7% 


321.75/332 
64.7% 



Table 1 : The first column lists the subset named as defined in the text together with its number 
of measurements and condition(s) if any. Each row displays the corresponding contribu- 
tion under the condition quoted in the title of the data column. The last row gives the total 
X^/(number of degree of freedom), followed by the fit probability. Information written bold- 
face indicates the distance of the fit function to a data set left outside from the fit procedure 
together with its number of data points. In this case, the condition parameter associated with 
the corresponding data set (scale or mass shift) is fixed to the value returned by the full global 
fit reported in the second data column and given in Tabled 

are always the improved errors returned by the routine MINOS. This has allowed us to check 
that the minimum was always locally parabolic, which provides symmetric errors. 

Tables [2] and [3] display the fit parameter values as returned by the fit under the conditions 
defined by the various titles of the data columns. We only provide results including the photon 
VP inside the definition of the pion form factor for e^e^ annihilations. Comparing the vari- 
ous data columns in Tables [2] and [3l clearly illustrates that the fit parameter values stay close 
together, and generally widely within their (MINOS) errors. The single exception is obtained 
while removing the spacelike data ; in this case, the coefficients for the subtraction polyno- 
mial li^J"' differ significantly from all other cases. This could have influenced the predicted 
values for the V — » e^e^ partial widths ; however, we have checked that this is not the case 
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Parameter 


Full Fit 


No r 


No Spacelike 


No p mass shift 


Scale New Timelike 


1.006 ±0.004 


1.000 ±0.003 


1.004 ±0.003 


1.007 ±0.003 


Scale Old Timelike 


1.012 ±0.009 


1.010 ±0.009 


1.011 ±0.009 


1.013 ±0.009 


Scale NA7 


1.008 ±0.007 


1.008 ±0.007 


1.008 


1.011 ±0.006 


Scale Fermilab 


1.006 ±0.007 


1.006 ±0.008 


1.006 


1.008 ±0.007 


CLEO Shift (MeV) 


0.40 ±0.52 


0.40 


0.36 ±0.52 


1.37 ±0.39 


6m^ (102 GeV2) 


-0.268 ±0.095 


-0.268 


-0.285 ±0.096 





a 


2.303 ±0.012 


2.297 ±0.012 


2.292 ±0.012 


2.306 ±0.012 


g 


5.576 ±0.015 


5.578 ±0.017 


5.597 ±0.016 


5.573 ±0.015 


X 


0.903 ±0.013 


0.902 ±0.013 


0.902 ±0.013 


0.903 ±0.013 


za 


1.503 ±0.010 


1.505 ±0.010 


1.507 ±0.010 


1.503 ±0.010 


zv 


1.459 ±0.014 


1.466 ±0.014 


1.453 ±0.014 


1.460 ±0.014 


zt 


1.246 ± 0.049 


1.245 ± 0.049 


1.243 ±0.049 


1.246 ± 0.049 



Table 2: Parameter values in fits performed including photon VP. Three data columns are 
associated with all data (first data column), removing only the r data (second data column) and 
removing only the spacelike data (third data column). The last data column reports parameter 
values returned while fitting all data sets by fixing 6m? = 0. Information written boldface 
displays values not allowed to vary in the fit procedure. 



numerically. 

Among the fit parameter values given in Table [2l the most interesting are clearly the fit 
values for the scale factors which are nicely consistent with the corresponding experimental 
information recalled at the beginning of this Section. 

The single exception is the CLEO global invariant mass shift which is found consistent 
with zero. Taking into account the way the 0.9 MeV expected shift has been determinecj^ [84], 
this information is interesting. It will be rediscussed when examining the fit residuals which 
provides an additional important information. 

The third and fourth data columns in Table [T] provide further information : 
i/ Removing the r data from the fit sample, one can construct the predicted distributions 
for the ALEPH and CLEO data which are fully derived using our model together with only 

^'^ In order to make consistent the p parameters derived from fits to the ALEPH and CLEO data. 



29 



light meson partial width decays and the e+e^ data. The predicted distance to CLEO data 
is practically unchanged with respect to fitting with them, while the prediction for the ALEPH 
data is not as good even if it remains reasonablelEl This indicates that CLEO data are in so nice 
agreement with predictions (especially e'^e~' data) that they do not really constrain the fit ! In 
this respect, ALEPH data, while introduced in the fit data set, clearly influence the procedure. 

ill Removing only the spacelike data looks a little bit more appealing. The distance of 
the NA7 and Fermilab data altogether is degraded by ~ 24 units, and the fit probability grows 
from 74 % to 93%, pointing to some slight difficulty in accomodating these data sets. However, 
this result is by no way problematic enough to either force us to remove the spacelike data or 
to deeply question their quality. 

In order to compute the distance of the data samples left out from fit, one had to choose 
either the CLEO energy shift or the NA7 and Fermilab scale factors, as they can no longer be 
fit. We choose to fix them to their fit value as given in the first data column in Tabled 

Tables [21 and [3l mostly illustrate that, whatever the fit conditions examined, the location of 
the minimum in the fit parameter space remains practically unchangecj^. The results obtained 
by removing r data from the fit sample, those by removing the spacelike data, those corre- 
sponding to removing the photon VP function (not shown) . . . are consistent with each other. 
Let us note that the fit parameter ci in the Ii^^{s) function {i.e. essentially the p self-mass 
function) has been fixed to zero, as it was not found to sensitively affect the fits in the energy 
range we are fitting. The various fit conditions only affect the fit quality which varies from 
good to very good, while including the photon VP. 

Figure [21 shows the fit with the e^e~ data in the timelike region superimposed. The 
global scale factor effects are accounted for. In the mass region, the lineshape is a predic- 
tion essentially derived from the phase and branching fraction of the (p tt+tt^ decay mode 
as measured by the SND Collaboration [71 J. Information on the full (local) invariant mass 
spectrum (when available) would certainly improve this prediction. 

Likewise, Figure [H shows the fit function and the ALEPH [[311 and CLEO |[33l data super- 
imposed. One may note that the highest data point from CLEO data lies at ~ 2a of the fitting 
curve. Actually, the lineshape of the CLEO data in the neighborhood of the maximum raises 
some difficulty while fitting, as will be seen shortly with the fit residuals. 

Leftside Figure [4l shows the spacelike data [[59l [80ll together with the fit function. One 
may note a small, but systematic upwards shift of the fit compared with the NA7 data which 
certainly explains the jump in the fit probability while removing this data sample from the fit. 

One may conclude from Figures [2l [3l [H and the fit probabilities that the agreement of the 
data with the model functions is good and that no obvious drawback shows up. 

In order to refine this statement, we have had a closer look at the fit residual plots. The 
fit residuals for the new Novosibirsk data are shown in Figure [51 The plotted experimental 
errors do not take into account the bin-to-bin correlations generated by the global scale error 
common to all data sets. As the difference between our pion form factor functions F^{s) and 

^^The distance for the fitted part of the ALEPH spectrum corresponds to an average per point of 1 .28 and 
thus to an average distance of l.lZa per data point. As a prediction, it is already a good starting point, which is 
improved by the global fit to an average distance of 0.90a per data point. 

^^With the exception mentioned above. 
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Parameter 


Full Fit 


No r 


No Spacelike 


No p mass shift 


Subtraction Polynomial : n^^(s) 


Cl 














C2 


-0.467 ±0.013 


-0.463 ±0.014 


-0.472 ±0.013 


-0.470 ±0.013 


Subtraction Polynomial : €2(5) 


Cl 


-0.071 ±0.003 


-0.071 ±0.003 


-0.072 ±0.003 


-0.072 ±0.003 


C2 


0.045 ± 0.004 


0.045 ± 0.004 


0.046 ± 0.004 


0.045 ± 0.004 


Subtraction Polynomial : ei(s) 


Cl 


-0.017 ±0.001 


-0.017 ±0.001 


-0.017 ±0.001 


-0.017 ±0.001 


C2 


0.020 ±0.001 


0.020 ±0.001 


0.020 ±0.001 


0.020 ±0.001 


Subtraction Polynomial n);^/'^(s) 


Cl 


0.918 ±0.061 


0.944 ±0.068 


0.727 ±0.074 


0.915 ±0.060 


C2 


0.433 ±0.106 


0.361 ±0.115 


0.831 ±0.145 


0.440 ±0.105 



Table 3: Parameter values under various strategies (cont'd). Boldface parameters are not 
allowed to vary. Each subtraction polynomial is supposed to be written ciS ± C2S^. 

F^{s) essentially lies in the p^ — uj — cj) mixing scheme produced by breaking Isospin Symmetry, 
this Figure can be compared with Figure 1 in ll26ll or Figure 9 in ll36l where a systematic s- 
dependent effect was pointing towards a consistency problem between e+e^ and r data ; such 
an effect is no longer observed, pointing to a presently more adequate manner of performing 
the breaking of Isospin Symmetry. One may possibly note that the dispersion of the residuals is 
very small everywhere for the 1998 CMD2 data (MIHII, while it is larger for the 1995 CMD2 
II27I and SND ll30l data which additionally, are moved in opposite directions in the uj mass 
region. This indicates that our fit parameter values are dominated by the 1998 CMD2 data. 

The residual distributions for r data -the upper plots in Figure [6]- look more interesting. 
The arrows indicate the limit of the fit regions. The errors plotted are certainly underestimated, 
as the bin-to-bin correlations are not accounted for in the drawings ; moreover, the errors 
produced by identifying invariant mass coordinate and central bin value also are not considered. 

One observes now a small but clear s-dependent structure above the p peak location (more 
precisely above s ~ 850 MeV) which certainly reflects the influence of the unaccounted for 
higher mass vector mesons. One has also examined the effect of removing the parameter 5m? 



31 



by fixing it identically to zero while fitting. The r data residuals are given by the lower plots 
in Figure |6l One clearly observes the rise of a structure at the p peak location in ALEPH data 
which is therefore a clear signal of a p° — mass difference. The mass shift observed is : 

m^o-mji = -5m'^ = (0.27±0.10) 10"^ GeV^ ^ m^o-mp± = 1.73±0.60 MeV (48) 

in good agreement with several other reported values ll54l . This question will be discussed with 
more details in the Subsection devoted the p meson parameters. 

In contrast, the shape of the CLEO residual distribution rather indicates a systematic effect 
in CLEO data located only in the p peak region. The global CLEO energy shift of 0.9 MeV 
serves to recover from the disagreement with ALEPH data. However, these plots clearly show 
that the problem of systematics is not global but local and that there is no evidence for a 
significant global invariant mass shift within the CLEO data in our fitted range. As this residual 
behaviour is also observed in the standalone fit performed in [f33l (see Figure 10 therein), it 
should not follow from the constraints specific to our model. 

Actually, there is a correlation between 5m'^ and the CLEO energy shift e which vanishes 
when performing a simultaneous fit of ALEPH and CLEO data. In order to check this state- 
ment, we have removed ALEPH data from the fit data set. Then, fixing 5m^ = 0, we get 
e = 1.57± 0.40 MeV and, conversely, fixing e = results in - 6m^ = (0.36 ± 0.08) 10^^ 
GeV^, with no change in the value and always the same residual shape as shown in Figure 
[6l Therefore, the value for 5m^ is set by the ALEPH data and should be confirmed by forth- 
coming data sets. Correspondingly, it is the use of ALEPH data which indicates that the CLEO 
energy shift could well be consistent with zero. 

All the reported exercises also show, as clear from Tables [21 and [3l that the fit parame- 
ter values are stable (with the exception already mentioned). This means that our model is 
overconstrained and that, practically, only the fit quality (i.e. the height of the minimum) is 
affected by the various conditions we have imagined. We were also aware of possible correla- 
tions between the subtraction polynomials. Looking at the fit error covariance matrix, we did 
not observe strong correlations between parameters belonging to different polynomials, which 
seems to indicate that they are indeed independent. 

Therefore, one may consider that the description of all form factor data supports our mixing 
model, as reflected by the statistical fit qualities reported in Table[i|under various conditions. 

Finally, rightside Figure |4] shows the predicted phase of the 1=1 part of the pion form 
factor together with the measured Pu (vrvr) phase data from If82ll83ll . Clearly, the description is 
good, keeping in mind that some contributions have not been included, especially the exchange 
of a spacelike p. Therefore, this Figure indicates that the neglected diagrams should contribute 
not more than a few degrees to the phase. 

Therefore, the description of all form factor data can be considered as satisfactory and 
provides a solid ground to our main assumptions : 

j/ The bulk of Isospin Symmetry effects which create the difference between e^e^ and r 
form factor lineshapes is a p° — a; — mixing scheme of dynamical {i.e. s-dependent) structure. 

jj/ An appropriate subset of meson partial width decays and the e'^e^ form factor data 
mostly suffice to set up a p*^ — — mixing scheme able to derive the r spectrum with good 
precision. 

jjj/ The effects of higher mass vector mesons in the mass region below 1.0 GeV, even if 
somewhat visible on the upper wing of the p peak, are negligible. 
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jv/ The (observed) p° — mass shift is very small and of negligible effect. New r data 
may confirm its relevance, as this follows only from ALEPH data. 

In view of the residuals shown in Figure [5] and, even more clearly in Figure [61 the 
statement in jjj/ can be precisely commented upon. Because of their statistical accuracy, the 
Novosibirsk data do not exhibit any undoubtful effect of higher mass resonances within the 
whole fitted range ; one may possibly guess a dip (of small significance) in the region between 
850 and 900 MeV. In the CLEO data (right plots in Figure (6]), the residual structure is unclear 
in the region between ~ 700 and ~ 950 MeV, while in the region above ~ 1 GeV, there is 
undoubtfuUy a missing structure which can reasonably be attributed to the p(1450) meson low 
mass tail. If one relies on the ALEPH data (left plots in Figure (6]), one is instead tempted to 
state that higher mass vector meson effects have influence down to ~ 850 MeV, {i.e. deep 
inside the high mass wing of the p(770) meson. Additional r data sets and large statistics e+e~ 
data sets (collected through the radiative return method) are expected in a near future from B 
factories and from DAPHNE, hopefully with controlled systematics. One may imagine that 
they should allow to clarify the situation in the region from the uj to peak. 

12.3 Light Meson Decays 

As a preliminary remark, when fitting partial widths (actually coupling constants), the rec- 
ommended data used are the partial widths taken from the RPP [54], when available. If not, 
they are derived from the branching ratios and the full widths. Sometimes, this procedure re- 
veals a surprising information. For instance, for r] — > 77, the ratio of the "fit" partial width 
error to the corresponding central value is 0.026/0.51 = 0.05, while the corresponding infor- 
mation derived from the quoted "fit" branching fraction is 0.26/39.98 = 0.007, which might 
look somewhat optimistic. 

The numerical estimates of branching fractions have been calculated using the information 
returned by the MINOS program and take into account the parameter error covariance matrix in 
the standard way (as recalled in Section 7.3 of [[T6l , for instance). This is mandatory as some 
error correlation coefficients are very large, namely those among the two fit parameters hidden 
inside ei(s), or inside £2(5) are about 95%. Most other error correlation matrix elements are 
below the 10 % level. We therefore consider that our error estimates are accurate. 

On the other hand, all partial width results we compute have been derived using the accepted 
values for all vector and pseudoscalar meson masses [|54l . In order to produce the branching 
ratios as given in Table |4l we have also divided these partial widths by the accepted total 
widths reported in the latest issue of the Review of Particle Properties [54J. The errors on 
masses and widths have been taken into account in the computer code used in order to derive 
the reconstructed branching ratios. 

12.3.1 Radiative Decays Of Light Mesons 

The fit values for the branching fractions of light mesons radiative decays are displayed on 
top of Table IH Most decay modes involving vector mesons are in nice correspondence with 
their recommended values [|54l . 
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The value returned for the u — > 7r°7 branching ratio is about 3 a from the presently rec- 
ommended value |[54l . but is in good agreement with the former recommended value ((8.5 ± 
0.5) 10"^ as well as the value found in the fit of the e+e" 7r°7 reported in [[7411 ((8.39 ± 
0.25) 10"^). This indicates that the recommended central value for this decay mode can be 
questioned. 

On the other hand, as already commented upon, the branching fraction found for u 777, is 
in much better agreement with the average value proposed by the PDG in [(54l| ((6.3±1.3) 10""^) 
than their so-called fit value reported in Table |4] ; this result is also in perfect agreement with 
the Crystal Barrel [|87l measurement ((6.6 ± 1.7) 10"'^) as well as the measurement obtained 
in a fit to the e+e~ r]'j cross section [|74l ((6.56 ± 2.5) 10~^). We also consider confidently 
our result for this decay mode. 

The ratio : 

T{uj rn) ^ 2 ± Q_Qll^ iQ-2 (49) 
T{uj 7r^7) 

depends only weakly on the mass and width definitions of the cu meson and is in agreement with 
all reported direct measurements in the RPP [.54il . This also gives support to both fit results. 

The contribution of the 77' p*^7 mode to the is 2.23, while all others are smaller or 
of order 1 . This may indicate that the box anomaly |[85l [88]| shows up and might have been 
accounted for. 

The only difficult point of the model is the ^ 1.9a departure of the partial width for t] — > 
77 from the expected value commented upon at the beginning of this Section. Whether this 
could be due to our assuming that the pion decay constant is not affected by Isospin Symmetry 
breaking effects is an open possibility. Instead, the partial width for r/ 77 fits nicely its 
expected value, possibly because of its larger experimental uncertaintjo 

Finally, we should note that our model gives a precise indirect measurement of fx/ fw '■ 



K 



1.503 ± 0.010,t,t ± 0.002wei ^ ^ = 1-226 ± O.OOAstat ± 0.001^^^^/ (50) 

Jtt 



where the second quoted uncertainty reflects details of the model together with the effects of 
including the spacelike data in the fit. This is in balance with the corresponding quantity which 
can be derived from the reported world average data [l54l as fx/ fw = 1.223 ± 0.010, assuming 
that the errors on f^, and /k are uncorrected. 



12.3.2 Leptonic Decays Of Light Vector Mesons 

Table m indicates that our model nicely accomodates the to e+e^ and e+e^ partial 
widths giving values which coincide with their recommended values Il54ll . 

Our result for p° e^e~ is derived from the same data which underly the other proposed 
values [54] and has been obtained with a careful account of all statistical and systematic re- 
ported errors. Therefore, this value can be confidently considered ; one should note that it 

In order to test this assumption, we have left /^o free in our fits and, for f^a = 87.9 ± 2.4 MeV, we have 
reached a probabiUty slightly above 80 %, with r(?7 — > 77) at 0.55 cr from its recommended value ||54| and 
T{jf 77) at only 0.37 a. This has to be compared with the reported value extracted from r(7r° — > 77) which 
provides /^o 91.92 ± 3.54 MeV. 
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exhibits a ~ lOa distance to the presently accepted branching fraction |!54l. A more straight- 
forward information coming out from our fits is the corresponding partial width : 

r(p e+e") = (8.34 ± 0.10 ± 0.31)10"^ MeV (51) 

where the first error merges all statistical and systematic uncertainties commented upon in the 
body of the text ; the second error takes into account the real uncertainty affecting the p mass 
used in order to derive the partial width from the coupling. It has conservatively been fixed 
to 10 MeV for reasons which will become clearer shortly. The corresponding partial width as 
given in L54J is (7.02 ± 0.11)10^3 ^eV. 



12.3.3 The uj/cj) ^ tt+tt Decays 

The value found for the (p tt+tt^ partial width compares well with its measured value 
ITtTI . Actually, one may suspect that this datum prominently influences some of our free pa- 
rameters, certainly those in the expression for ei(s). The phase of the corresponding coupling 
constant being close enough (1.4(t) to expectation [71 J might indicate that the data (modulus 
and phase) for this mode carry small systematic uncertainties. 

The branching fraction we get for the uj tt+tt^ mode is more appealing. It is derived 
from all data involved in this measurement with a precise account of all systematic uncertain- 
ties. Additionally, the quality of the measurement we propose probably does not suffer from 
significant model uncertainties, as the p — uj interference region is quite well described (see the 
insets in Figure O. Therefore, our conclusion for this decay mode is either of : 

Br(c<; ^ TT+TT-) = (1.13 ± 0.08)% , V{uj tt+tt-) = (9.59 ± 0.80) 10^^ MeV (52) 



using the recommended value for width and the mass of the u meson 11541 . 

This new datum may influence the global fit of all the cu decay modes in isolation. This is of 
concern for our purpose, as one has noticed that the disagreements observed between the PDG 
recommended values [54J and our results for the VP^ modes refer mostly to the uj ^ [r]/ 7r°)7 
branching ratios. Along this line, our fit solution provides : 

^ = 0.14 ±0.01 (53) 

1 (cu — > 7r^7) 

to be compared with the single existing measurement by KLOE fSQ^ which provided 0.20±0.04 
and a 1.5(T distance. Therefore, our surprising estimate of the two pion mode together with 
radiative decays is in accord with experimental expections. 

Finally, the unfitted Orsay phase for the coupling uj tt+tt" is found close to its expected 
value from a standalone fit to the so-called old timelike data [50|, while our fit for the phase of 
the TT+TT^ coupling is in good agreement with its measured value fT\\ . 



12.3.4 The KK Decays 

As explained in the body of the text, we have been led to leave both — > KK de- 
cay widths outside our fit procedure, as there is some uncertainty with possible factors, like 
Coulomb corrections, which may affect the usual coupling constant contributions to both par- 
tial widths flTl |73l |90l . 
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Therefore, the values reported in Table |4] are predictions only influenced by the other de- 
cay modes and without any additional correction factor to each of the — > KK branching 
ratios. The numerical values found for these branching ratios clearly illustrate that our model 
is overconstrained and provides precise values for the coupling constants of both KK 
modes. 

The ratio of the prediction to the recommended central values is 1.022 for the charged 
mode {i.e. a 1.8 a distance) and 0.97 for the neutral decay mode (i.e. a 2.0 a distance). Taking 
into account the model uncertainty reported in Table IH the agreement could be considered as 
satisfactory. 

Now, if correction factors have to be applied, they are expected to improve the prediction 
for the rates. Therefore, they should be of the order 0.976 and 1.031 for, respectively, the 
charged and neutral decay widths. This clearly invalidates the traditional 1.042 correction 
factor proposed in order to account for Coulomb interaction among the charged kaonQ 

Correcting both modes as argued in 11731 . even if able to provide a good account for the ratio 
r(0 K^K^)/r{(f) — > K'^K^), does not allow a good account of both modes separately, as 
the corrections proposed turn out to increase the expected rate for both modes. 

Within the framework of our model, if correction factors have to be applied, they should not 
increase the charged decay mode by more than ~ 1%. There is more freedom with the neutral 
decay mode. Therefore, in order to fix one's ideas, one has let a correction factor for only the 
neutral decay mode to vary. In this case, of course, the correction factor (= 1.047 ± 0.024) 
is found such that the neutral mode exactly coincides with its measured value, which could 
be expected beforehand. However, more interesting is that the contribution of the charged 
decay mode (which does not explicitely depend on this factor) is only 0.3 (a 0.5 a effect). 
This indeed confirms that only the predicted neutral decay width might have to be corrected 
significantly. Taking into account that systematic effects are harder to estimate for the charged 
mode than for the neutral one (see footnote 16 in [16j|), this may look a physical effect. Whether 
the "mixed isoscalar and isovector source" scheme of f90] can account for such an effect would 
be interesting to explore. 

As a summary, our analysis tends to disfavor a significant correction factor to the — > 
K^K^ decay width (above the 1.01 level). It would rather favor a significant correction factor 
for only the neutral mode K^T(^ (which could be as large as 1.047 for the rate). If 
the traditional scheme of Coulomb corrections should really apply, both measured widths for 
— > KK are hard to understand, as already stated in 11721 [T6ll . 

12.3.5 What Are The p Parameters? 

For objects as broad as the p (or K*) meson, the definition for mass and width (from ex- 
perimentally accessible information) is not a trivial matter [|53l , and no physically motivated 
uniqueness statement can be formulated. Having defined in our model the p° and p"^ propa- 
gators as analytic functions (or rather meromorphic functions on a 2-sheeted Riemann surface 
with branch point at threshold), one has at disposal the poles of the propagators. This has been 
shown to provide the most stable definition of the mass and width ||53l . If one assumes that 
Analyticity of S-matrix elements is a basic principle, this is also the most model independent 

See ll72l for a detailed account of the usual way to deal with Coulomb corrections and Isospin Symmetry 
breaking effects in (p decays. 
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definition. Indeed, wliatever are tlie working assumptions, the pole basically tells where the 
peak is really located and how wide is the invariant mass distribution around the peak (typi- 
cally, close to the full width at half maximum). Obviously, any given model cannot successfully 
describe the relevant data without having the pole located at the place requested by the data. 
The specific character of a given model is basically concentrated in the regular part of the Lau- 
rent expansion of the amplitude. To be complete, departures from this statement may exist but 
is (parameter, model) freedom is essentially governed by the magnitude of the experimental 
error bars. 

However, based on the expectation that the mass difference between charged and neutral p 
mesons is only of the order a few MeV's at most {i.e. Snip ~ a few 10~^mp), one may guess 
that this difference could be less sensitive to mass definitions. 

Our final results for the complex s locations of the p meson poles can be derived from our 
fit parameter values. Sampling them by taking into account the parameter error covariance 
matrix, one gets : 

' Spo = (0.5782 - i 0.1099) ± (0.9 + i 0.5) lO^^ (GeV^) 
< Sp± = (0.5760 -i 0.1095) ±(1.0 + i 0.5) 10-=^ (GeV^) (54) 

^ - Sp± = (2.26 - i 0.38) 10"=^ ± (0.83 + i 0.14) 10"=^ (GeV^) 

with uncertainties folding all reported statistical and systematic errors. 

In order to compare with related information available in the literature, one has to relate the 
p pole locations with the usual Mp and Tp. 

Defining as ll56l sr = M| — iMrTr, one gets : 

' Mpo = 760.4 ± 0.6 MeV , T^o = 144.6 ± 0.6 MeV 

< Mp± = 758.9 ± 0.6 MeV , Tp± = 144.3 ± 0.5 MeV (55) 

^ MpO - Mp± = 1.51 ± 0.53 MeV , T^o - Tp± = 0.22 ± 0.08 MeV 

which for p° are slightly larger than those found by ||56l using only the so-called old timelike 
data. One may note that the mass difference is affected by a smaller uncertainty than masses 
separately ; this effect is even much more pronounced for the width difference. This error 
shrinking reflects the correlations contained in the parameter error covariance matrix of our fit. 

One may also choose ll57ll sr = {Mr — ir/j/2)^ and obtain slightly different values (not 
provided). Defining the mass by the location of the maximum of the distribution and the width 
by the full width at half maximum cannot be derived easily from Eqs. (|54l ) ; they are : 

' MpO = 762.1 ± 0.6 MeV , TpO = 144.5 ± 0.6 MeV 

< Mp± = 760.8 ± 0.6 MeV , Tp± = 144.5 ± 0.5 MeV (56) 

MpO -Mp± = 1.22 ±0.53 MeV , r^o - T^i = 0.02 ± 0.08 

One may consider these values as they are in consistency with the way the mass and width for 
objects like the cu and mesons are usually defined ||54l . The difference between these results 
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and those in Eqs. (l55l) could be attributed to the influence of the regular part of the invariant 
mass distribution which distorts a little bit the distribution lineshape. Finally, in view of Eqs. 
(l55l) and (l56l) . one cannot be really conclusive about the sign of the width difference central 
value, as it sensitively varies with parameter definitions (see also footnote (6]). 

Our mass difference values can be compared with results from other reactions available 
in the literature [.54,1 . Limiting oneself to the most recent estimates (all with large statistics). 



one has : 












nipQ 


- mp± 


= 0.4 ±0.7 ±0.6 


MeV 


KLOE 2003 [89], 1980 Kevents 


< 


TUpO 


- mp± 


= 1.3 ±1.1 ±2.0 


MeV 


SND 2002 m], 500 Kevents 






- mp± 


= 1.6 ±0.6 ± 1.7 


MeV 


Crystal Barrel 1999 [92], 600 Kevents 



which compare satisfactorily with our fit results in either of Eqs. (|55l) or Eqs. (l56l) . These 
experiments analyze the 7r+7r~7r° final state produced in e+e~ annihilations at the mass II89[ 
[91 J or in annihilations at rest [92J using standard varying width Breit-Wigner shapes for 
both charged and neutral p meson distributions. 

The ALEPH Collaboration has also performed a global fit ||3TI of the e~^e~ Novosibirsk 
data together with the ALEPH and CLEO r data sets (as we did) and gets : 

rripo - mp± = -2.4 ± 0.8 MeV ALEPH 2005 [3l] (58) 

significantly different from our findings. However, their fit residuals show a s-dependence 
below s ^ 1 GeV^ which is absent from our datc|^. 

The fit results from other experiments shown in Eqs. (|57] ) are consistent with either sign for 
nipo — mp± ; our own results favor nipo — mp± > (with resp. 2.9 a and 2.3 a significance), 
while naively one may expect the opposite. However, Bijnens and Gosdzinsky ||93l , analyzing 
within the ChPT framework all contributions to this mass difference, concluded that : 

- 0.4 MeV < rripO - mp± < ±0.7 MeV (59) 

All measurements given in Eqs. (1571) are consistent with this mass interval. The ALEPH mass 
difference is at a 2.5 a distance from the lower bound of Eq. (|59l ). Concerning our results, our 
largest estimate of the mass difference (1.51 ±0.53) is only at 1.5 cr from the upper bound, while 
our smallest estimate (1.22 ± 0.53) is at 1.0 a from this upper boun^^. Therefore, awaiting 
new measurements which may confirm the ALEPH spectrum at the peak location, one may 
conclude that our fit results are in good agreement with ChPT expectations. From a statistical 
point of view, this agreement has even been marginally improved compared with the ALEPH 
best fit result. 



^''Ref. ||3T1 does not provide, strictly speaking, the fit residuals, but this can be guessed from Figure 67 therein. 
^"Comparing with the upper limit is the most favorable case for our results, while comparing with the lower 
limit is the most favorable case for ALEPH. 
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If one decides to parametrize the distributions with varying width Breit-Wigner shapes, 
one would recover more traditional mass values as tabulated in [|54l . However, their model 
dependence (not only their definition dependence) should be stressed. Because of having decay 
and pion form factor data intricated, it is not easy to perform this check within the present 
context. However, a good approximation of using varying width Breit-Wigner expressions is 
to define the p masses by : 

Re{D-^\Ml)) = (60) 

which, by the way, is fulfilled by the standard Gounaris-Sakurai propagator [|94l [311 . In this 
case we get : 

Mpo = 774.8 ± 0.6 MeV , M^± = 773.3 ± 0.7 MeV , MpO - Mp± = 1.48 ± 0.50 MeV 

(61) 

which can be compared with standard values for masses and provide a mass diffference in 
good agreement with our results above. This definition of p meson parameters is very close to 
the corresponding ones following from using Gounaris-Sakurai expressions. Then, comparing 
Eqs. (|6T1) with the corresponding ALEPH results Oil is a way to exhibit the effect of the 
additional isospin breaking mechanism we propose. 

Going a step further along the same line, one may also choose to define masses by the HK 
mass values as they come from the Lagrangiarl^l] and our fitted parameters. Then one gets : 

MpO = 782.1 ± 2.1 MeV , Mp± = 780.4 ± 2.2 MeV , MpO - Mp± = 1.74 ± 0.60 MeV 

(62) 

This result is also interesting. Indeed, as stated above (see Footnote (6]), a reasonable break- 
ing of Isospin Symmetry at Lagrangian level, while producing a (HK) mass difference between 
the p° and p^ mesons returns rupo = m^. As one may think that the HK cu mass is close to 
the tabulated value [54] (m^ = 782.65 ± 0.12 MeV), it may have a meaning to find that the 
HK mass for the p° is consistent with the accepted u mass. For this purpose, it should be noted 
that the u mass value used in our fits was fixed at this accepted value, and then cannot directly 
influence the HK value for rUpO . 

Defining the width using Eq. (II 141) and usindj the HK p masses given in Eqs. (|6T1) one 
derives the following information for the p widthO : 

Tpo = 174.3 ± 2.0 MeV , Tp± = 175.0 ± 2.0 MeV , T^o - Tp± = -0.7 ± 0.2 MeV. 

(63) 

which are substantially larger than all other definitions and with a difference going in the oppo- 
site way compared with above ; however this follows from the most usually employed formula 
for the two-body decay widths of the cu and </> mesons. One may argue that the larger mass and 
width exhibited by Eqs. ([621) and ([631) compared with usual is related with our having the HLS 
parameter a significantly different from 2 (a 3 cr effect). 

This has a clear physical meaning : It is the mass of the p meson while working at tree level (when possible, as 
maybe farther inside the spacelike region than we have gone). However, one should keep in mind that, numerically, 
they are derived from fits using expressions containing the unavoidable self-mass corrections. 
^^We have averaged the pion masses used for the charged p width computation. 

We tried to decouple gpO^^ from 5p±7r7r by defining gp±T,Tr = gp«7nT + ^9 with 6g submitted to fit. The 
minimization procedure returned 5g — —0.018 ± 0.014 and an improvement of the of the order 0.5, quite 
unsignificant. 
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In summary, as far as p meson parameters are concerned, we consider that the most relevant 
information are the s location of the poles as given by Eqs. (|54|) ; definitions of the mass and 
width using sr = — iMbTr, or sr = {Mr — iTR/2)'^, or something else, can be derived 
algebraically and some have been given. 



12.4 A Few More Comments On The Model 

In order to justify the change from ideal to physical fields, one should check that the func- 
tions in the non-diagonal elements of the R matrix in Eq. (fTTI) are small compared to 1 in the 
relevant invariant mass range. For this purpose, we have had a closer look at the functions : 



n^Trfs) - £2(5) 



Fpd^)= n , 2 fn , ^ TTT ^64) 



[1 - zv)m? + (1 - /i^)e2(s) 



computed with the fit parameter values. These are the entries of the R{s) matrix which defines 
our transformation from ideal (bare) fields to physical fields. 

Figure |7] shows the real and imaginary parts of these functions. They are all small compared 
to 1 all along the physical region : \FpJ{s)\ ~ 0{lQ-'^), \Fp^{s)\ ^ C(10-2) and \F^^{s)\ ~ 
(9(10^^). As expected F(^(/,(s), which does not vanish in the Isospin Symmetry limit, is larger 
and one observes an order of magnitude difference. 

Fpu{s) represents the traditional p — uj mixing and its behaviour translates in our modelling 
the known large (Orsay) phase by the quasi-vanishing of its real part around m^. It exhibits 
around ^ ^ 0.3 GeV, the two-pion threshold, an unexpected behaviour which is actually 
too small to influence numerically the pion form factor. Whether this local effect should be 
considered seriously is unclear, taking into account the approximations done in order to work 
out the mode0. The important point here is that, even if narrow, the amplitude does not exceed 
a 10 % level there and a few percents all along the physical region (namely the p/uj and (p peak 
locations). 

Fp^(s) is of the same order of magnitude than Fp^(s) but much smoother all along the 
physical region ; its real and imaginary parts around the (j) meson mass are comparable (~ 1 
%). 

Fu)(f){s) is more interesting as it represents what is traditionally attributed to a constant 
mixing angle of a few degrees [[SJ. It is indeed what is exhibited, as this function is real at 
and close to real at m^. However, the numerical values of the mixing angle vary significantly : 
at the tu mass the angle is 0.45°, while it is 4.64° at the mass. This tends to indicate that the 
notion of mixing angle has somewhat to be readdressed. 



^'^This goes along with the remark that one would have preferred a solution for Fp^{s) which vanishes at s = 
as Fpci,{s) and Fa,0(s) do, even for non identically vanishing of ei. This could also be a consequence of working 
at first order in ei and £2- However, the neglecting of the anomalous loop effects may play some role near s — 0. 
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As explained above, the (not too small) magnitude of F^^{s) could have been inferred from 
the HLS model, as the transitions between luj and 0/ follow from (kaon, K*K and K*K*) loop 
effects and not from supplementing them with Isospin Symmetry breaking effects. 



Figure [8] emphasizes the important features of our model for the pion form factor. The 
upper plot shows the function : 

which summarizes the breaking of Isospin Symmetry all along the physical region. The strong 
effects at the to and mass locations could have been expected. However, one clearly sees a 
non-zero "background" contribution extending down to threshold (and even below) and beyond 
the (p mass. This simply illustrates that our p — cu — (p mixing scheme is really invariant mass 
dependent. This is why it can suppress the unwanted effects exhibited by residuals in more 
standard approaches (see for instance [|26ll , [|35l or If36l ). 
The lower plot shows instead : 

2 II?r/„M2 



s 



Hi{s)= ' ----- (66) 

where F^ j^^(s) is identified with the p° part of the pion form factor, as traditionally done. It 
clearly exhibits that the p° and mesons are different kinds of objects in our modelling. We 
indeed observeEfl an effect of several percents and functionally s-dependent. 

Stated otherwise, the is indeed a pure isospin 1 meson, while the p° meson is actually 
a mixture of isospin 1 (pj) and isospin states (cuj and (pj). A real extraction of the isospin 
1 component of F^{s) should isolate the isospin 1 part of the p, cu, and </> amplitudes, which 
is exactly what our model does, by construction. This allows us to agree with the analysis by 
K. Maltman Ii48ll49l concluding that the p part of the pion form factor in e+e^ data does not 
behave as being isospin 1 ; however, this does not invalidate the e+e^ data. 

Therefore, the single departure from CVC one observes is simply the tiny shift between the 
p*^ and masses which follows from ALEPH data. 

It also follows from our model that using r data in order to reconstruct the equivalent e"'"e~ 
spectrum is not straightforward and requires a non-trivial physics input as shown by both plots 
in Figure [8l 

One is tempted to think to extending the model in order to include higher mass vector 
meson nonets. In order to do it properly within the HLS context, one may change in the co- 
variant derivative (see Eqs. (1731) ) the (presently) single vector term gV to 91^1- This induces 
further transition amplitudes like, for instance, p°(770) < — > p°(1450) which may sharply com- 
plicate the model. It is not easy to have a feeling for the numerical magnitude of the inter-nonet 
transitions as, now, pion and kaon loops contribute. Moreover, the higher mass vector mesons 
are in the region where all the thresholds of the VP channels are open. Therefore, one may 



^^The local minimum just above threshold, which reflects the structure of Fpui{s) we discussed, may be an 
artifact of the model. Its numerical value at this location makes it however totally invisible. 
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also have to include the corresponding loops explicitly. Whether the problem will not become 
numerically untractable is therefore unclear. 

For this pupose, one has certainly noticed that the problem we have examined is highly 
non-linear in the parameters. This makes the search for solutions highly dependent on getting 
a starting point (in the parameter phase-space) reasonably close to the final solution that the 
fit procedure may succeed. In the present case, it was already a non-trivial (and highly time 
consuming) task. 

Finally, as clear from the downmost Figure [H our model predicts a few percent effect for 
isospin breaking in the resonance region as well as in the low energy region where ChPT 
applies. Although our pion form factor fit results are good down to the lowest energy points 
in e+e^ data (about 360 MeV), it would be interesting to compare with ChPT estimates in this 
region. Examining this question, the authors of Ref. Il23l argued that a reliable answer may 
call for a 0{p^) estimate still missing. 

13 Conclusion 

Within the context of the HLS model, as for other models, at tree level the so-called ideal 
fields (p/, toj, (f)j) are mass eigenstates. This simple picture vanishes at one-loop order. In this 
case, kaon loops generate non-zero amplitudes allowing tu/ < — > 0/ transitions. Breaking of 
Isospin Symmetry in the pseudoscalar sector generates a mass difference between kaons and, 
besides, the transition amplitudes cui < — > pi and (pi < — > p/, even if of small magnitude, are no 
longer vanishing. Additionally, these amplitudes exhibit a dependence on s, the square of the 
momentum flowing through the vector meson lines. As the physical fields p, cu, (p are expected 
to be eigenstates of the squared mass matrix, this unavoidably leads to define them as linear 
combinations of their ideal partners. However, as the transition amplitudes are s-dependent, it 
is clear that these combinations should also be s-dependent. 

We substantiate these considerings starting from the HLS Lagrangian, modified by includ- 
ing in the squared mass matrix of the neutral vector mesons all self-energies and transition 
amplitudes. Making the assumption that the physical neutral vector fields should be eigen- 
states of the loop modified squared mass matrix of the (ideal) neutral vector meson, we solve 
the eigenvalue problem perturbatively. This leads to physical vector meson fields expressed 
as linear combinations of their ideal partners with definite s-dependent coefficients, which are 
actually analytic - or, rather, meromorphic - functions of s. Of course, this algebra leaves 
unchanged the charged and/or open strangeness sector, as the starting fields acquire a running 
mass but no transition from one to another meson field. 

The main mechanism producing the vector meson field mixing is the occurence of neu- 
tral and charged kaon loops in transitions between the ideal neutral vector meson fields. We 
have also shown that the anomalous HLS-FKTUY sector and the Yang-Mills lagrangian piece 
provide as supplementing mechanisms the K*K and K*K loops, occuring in transitions in 
nice correspondence with kaon loops. Within this framework, the cucp mixing has been shown 
to proceed from quantum (loop) effects, while the p^to and p°0 mixings follow from Isospin 
Symmetry breaking effects and vanish when this symmetry is restored. From a numerical point 
of view, the bulk of the effect is carried more by the subtraction polynomials than the loop 
expression themselves, as these are logarithmic functions of small amplitude in the physics 
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region of interest. 

The vector meson mixing allows to dynamically generate isospin violating couplings (pn'^n' 
and a;7r+7r~ at the (modified HLS) Lagrangian level. With this at hand, we have been able to 
construct the pion form factor expression at one-loop order modified in order to account for 
Isospin Symmetry breaking through only a dynamically generated — cjj — mixing scheme. 
A priori, this fully affects the e+e" physics but, by no way, the r physics ; in this sector, as 
a second order refinement, we have been led to accept shifting with respect to each other the 
Higgs-Kibble masses of the p° and mesons. This provides a tiny effect, nevertheless clearly 
visible in ALEPH r data, but not obvious in CLEO data. New data expected from B -factories 
may confirm the need for this mass shift. As a final result, we end up with structureless residu- 
als in the fitted regions, which confirm that our dynamical mixing scheme is appropriate. 

The mixing properties have been introduced in the anomalous decay amplitudes V — > P7 
and P — > V^. These processes actually represent our main lever arm while defining numeri- 
cally our isospin breaking scheme. 

Beside a good description of the V P7 and P decay data, this allows a very 

good simultaneous description of all pion form factor data from the close spacelike region to 
the mass, in e+e~ annihilations as well as in r decays. The physical ground of this result 
can be traced back to the fact that the p° meson is a mixture of isospin and 1 states (as the 
and uj meson), in contrast to the meson which is purely isospin 1. Actually, extracting the 
isospin 1 part of the pion form factor in e+e^ annihilations requires to split up the and cu 
contributions. This is done automatically by our model, and one can claim that such a splitting 
cannot be done without some model. 

The net result of this model is to prove that the lineshapes for the pion form factor in 
e~^e~ annihilations and in r decays are perfectly consistent with each other, without any further 
breaking of CVC than a possible tiny —p^ mass difference. In a further step, one may include 
the data from KLOE as well as data expected to come from BaBar and Belle concerning the 
r spectrum on the one hand and the e+e Initial State Radiation samples on the other hand. 
However, our comparison of e+e~ and r data does not seem to leave room for any kind of new 
physics. 

On the other hand, we have shown that this model allows a good account of all decays 
of the form V — > P7 and P V'y. The case of the a; — > 777 and u — > 7r°7 partial widths, 
where some disagreement is observed with the so-called "fit" values proposed by the PDG, has 
been discussed and we have argued that the real situation is somewhat unclear. 

Our dealing with the pion form factor data has led us to propose improved values for data 
sharply related with the e~^e~ tt+tt^ annihilation process, namely the p° e+e^ and 
tu — s> TT+TT" partial width decays. In both cases, we find that the reference values should be 
significantly modified and we propose for these new reference data. 

Finally, we have briefly commented upon the mass and width of the p° and p"^ mesons 
and argued that the best motivated definition should rely on the pole position in the complex 
s-plane and related definitions. 

Having shown that e+e" data and r data are perfectly consistent with each other pro- 
vided one uses an appropriate model of Isospin Symmetry breaking, we can conclude that there 
is no reason to question the e+e" data. This result is important as these data serve to estimate 
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numerically the hadronic photon vacuum polarization used in order to predict the value of the 
muon anomalous magnetic moment g — 2 . Therefore our model indirectly confirms the 33 a 
discrepancy between the BNL direct measurement of the muon anomalous magnetic moment 
and its theoretical estimate. 
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Decay Mode 


Fit Value 


PDG/Reference 


p [xlO"^! 

r / L J 


5.17 ± 0.04 


6.0 ± 0.8 


p 7r=^7 [xlO^l 

r / L J 


5.03 ± 0.03 


4.5 ± 0.5 


^ rry [xlO^l 

r 7 / L J 


3.05 ± 0.04 


2.95 ± 0.30 


ri' p-i [xlO^l 

/ III J 


33.3 ± 1.0 


29.4 ±0.9 


i^*± ^ i^'±7[xlO'^] 


9.8 ±0.9 


9.9 ±0.9 


1 I J 


2.26 ±0.02 


2.31 ±0.20 


^ ^ 7r°7 [xlO^] 


8.23 ± 0.04 


8.9+°-2^ (*) 


uj — > ?77 [x 10^] 


6.60 ± 0.09 


4.9 ± 0.5 (*) 


7]' io-y fx 10^1 

/ ' L J 


3.14 ±0.10 


3.03 ±0.31 


^ 7r°7 [xlO^] 


1.24 ± 0.07 


1.25 ± 0.07 


(h^rn fxlO^l 


1.292 ± 0.025 


1.301 ± 0.024 


^ r7'7 fx 10^1 

r III J 


0.60 ±0.02 


0.62 ±0.07 


r? — > fx 10^1 

'/ / / L J 


35.50 ± 0.56 


39.38 ± 0.26 


77' ^ 77 fx 10^1 

/ III J 


2.10 ±0.06 


2.12 ±0.14 


p e+e" fx 10^] 


5.56 ± 0.06 


4.70 ± 0.08 (**) 


^ e+e" fx 10^] 


7.15 ± 0.13 


7.18 ± 0.12 


0^ e+e-fxlO^l 


2.98 ±0.05 


2.97 ±0.04 


^ TT + TT^fxlO^] 


1.13 ± 0.08 


1.70 ± 0.27 (**) 


Qum+w- phase [degr] 


101.2 ± 1.6 


104.7 ± 4.1 (**) [|50l 


^ TT + TT^fxlO^] 


7.14 ± 1.7 


7.3 ± 1.3 


qstt+tt- phase [degr] 


-27.0 ±0.5 


-34 ±5 [Ell 


^ i^'+i^'-fxlO^] 


50.3 ± 1.0 


49.2 ± 0.6 (**) 


0^i^OirOfxlO2] 


33.0 ±0.7 


34.0 ±0.5 (**) 



Table 4: Reconstructed Branching fractions for radiative and leptonic decays using any of the 
various fit strategies. The reported values are the mean value and the rms of the simulated 
distributions. The last column displays the recommended branching ratios (5^. The symbol 
(*) indicates data commented upon in the text, (**) indicates data which are not introduced in 
the fit procedure. 
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Appendices 



A The Full HLS Non-Anomalous Lagrangian 

The construction of the non-anomalous Lagrangian of the Hidden Local Symmetry (HLS) 
Model has been presented in great detail several times by its authors (see for instance [UOl or 
more recently [2J). Let us simply outline the main steps of the construction procedure. 

The HLS model allows to produce a theory with vector mesons as gauge bosons of a hidden 
local symmetry. These acquire a mass because of the spontaneous breakdown of a global chiral 
symmetry Ggiobai = U{3)l^U (3)r. The chiral Lagrangian is written : 

>Cchirai = ^Trid^Ud'^U] (67) 

where U (x) = exp [2iP(x) / /tt] ; here is identified with the pion decay constant (/^ = 92.42 
MeV) and P is the matrix of pseudoscalar mesons (the Goldstone bosons associated with the 
spontaneous breakdown of GgiobaO- This matrix : 



V2 



7S* 



71"° 



TT 



v 



TT 



K- 



K 



(68) 

contains a singlet term besides the octet term ; appropriate combinations of T^g and t]q corre- 
spond to the physical pseudoscalar fields r] and r]'. Here and throughout this paper we restrict 
ourselves to three flavours. 

However, besides the global symmetry Ggiobai, the chiral Lagrangian possesses a local sym- 
metry iYiocai = SU{3)v which is included in the HLS approach by splitting up f/ as : 

U{x) = (69) 
where the ^ fields undergo the local transformation. These variables are parametrized as : 

^ia/fa^±iP/fn 



R,L 



(70) 



and the scalar field a is usually eliminated through a gauge choice, and can be considered ab- 
sorbed into the gauge bosons and removed. However, the decay constant f^- goes on appearing 
in the model through the HLS fundamental parameter a = fl/ fl- Using this parametrization 
Eq. (|67l) can be rewritten : 



chiral 



t ^21 



(71) 



This Lagrangian can be gauged for electromagnetism, weak interaction and the hidden local 
symmetry by changing the usual derivatives to covariant derivatives -D^ [[ll El and one then 
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gets : 



HLS 



La + aL.\ 



4 



(72) 



using obvious notations. 

Now let us turn to the covariant derivatives. These are given by[|2l : 

(where we have factored out the universal vector meson coupling constant) with : 



(73) 



L^ — cQA^ + 



cos 



^^(T, - sin^ ew)Z, + + W~T.) 



n. 



eQA, 



92 



(74) 



■sin 9]yZ^ 



cos 9w 

Eqs. (1731 1741) introduce the matrix of vector meson fields (the gauge bosons of the 
hidden local symmetry) which is : 



V 



V2 



P 



K*+ \ 



V 



p 

K* 



¥ ) 



(75) 



in terms of the so-called ideal field combinations (indicated by the superscript /) for the neutral 
vector mesons, which should be distinguished from the physical fields introduced in the main 
text. is the electromagnetic field and e the unit electric charge, g2 and 9w are respectively the 
gauge weak coupling constant and the weak (Weinberg) angle. and W^^ are, of course the 
weak gauge boson fields. Q, the quark charge matrix and the T matrices are SU(3) matrices : 



Q = 1/3 Diag(2, -1, -1) and 



-1, -1) , while T+ = (r_)t with : 




(76) 



1/2 Diag(l, 

/ Vud 

V 

in terms of elements of the Cabibbo-Kobayashi-Maskawa matrix elements. In this work, the 
terms have not to be considered ; they have been given for completeness. The HLS Lagrangian 
given above should be completed with the vector meson kinetic energy term[2] but also with 
the usual free Lagrangian for electromagnetic and weak boson fields. The leptonic sector also 
has to be added ; it is written as per usual : 



E 



e={e,tJ.,T) 



92 

2V2 



(77) 
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From a practical point of view, g2 defined above is related with the Fermi constant Gp and 
the W boson mass by : 

g2 = 2mw \/GfV2 (78) 
and it is useful to note that at the r lepton mass scale one has [54] : 

g2 = 0.629 (and e = 0.30286) . 

B The HLS Anomalous Sector 

QCD admits a non-abelian anomaly which explicitly breaks chiral symmetry. This anomaly 
is well reproduced by the Wess-Zumino-Witten Lagrangian f95^, ; this has been incorpo- 
rated within the HLS context by Fujiwara, Kugo, Terao, Uehara and Yamawaki along with 
vector mesons [I3l|2l. In this way, it becomes possible to provide a framework which allows 
one to describe most decays of vector mesons, and especially modes like u tt+tt^tt^ and 
others more relevant in the present context. 

Let us briefly outline the derivation and its assumptions which has been presented in com- 
prehensive reviews ||71|971|2l|. The anomalous action can be cast under the form : 

r = Fwzw + Tfktuy 

(79) 

Tfrtuy = J2i=i Q / d^x Ci 

where Fwzw is the original WZW Lagrangian. The Lagrangian pieces Ci where first given in 
|[3l and each of them contains APPP and AAP pieces which would contribute to the anoma- 
lous decays, but are cancelled by APV terms. These Lagrangians contain also VPPP and 
VVP pieces E |2l|- A priori, the weighting coefficients q are arbitrary. However, in order 
to reconcile this Lagrangian with decay data, especially uj tt+tt^tt^, FKTUY [l3l|2l finally 
choose the following combination : 

1 

^FKUTY ^ _^^,/..p.Tr[9^K9,KP] - -C,ppp , (80) 

which turns out to complement the usual WZW term for 'jPPP interaction with only a VVP 
term. In this model, for instance the decay tt^ — > 77 occurs solely through 7r° — > tup*^ followed 
by the (non-anomalous) transitions c<j 7 and p° — ^> 7 and the partial width is identical to the 
Current Algebra prediction reproduced by C^^p : 

TV p2 

C,,p = --^e^'P^'d^A^dpA^TilQ^P] (81) 

In 

The model given by Eq. (|80l) with : 

ie N 

C,ppp = --^e^''''''A,TT[Qd,PdpPd^P] (82) 
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(Nc = 3) has been also shown [85 J to reproduce perfectly the data on rj/rj' — > 7r+7r^7, es- 
pecially the most precise ones [88J. Indeed, with no free parameter, the distortion of the p 
lineshape is accurately accounted for and this can be considered as the signature for the box 
anomaly in experimental data. 

Accounting for the light meson radiative decays of the AVP or AAP forms is then a VVP 
coupling followed by one or two V — > 7 transition(s). From a practical point of view, it has 
been also shown [9J that the corresponding couplings can be directly derived from the following 
Lagrangian piece : 

C = Ce>^^p''i:i[d^{eQA, + gV,)dp{eQA, + gV,)P] , C = . (83) 

Let us note that in the meson decays we are interested in, the weak boson sector is irrelevant. 
Finally, one can find the VV P Lagrangian expanded in [f8l|, more precisely in Appendix 1 and 
4 for respectively the fully flavor symmetric case and the SU(3)/U(3) broken casCj 

C SU(3)/U(3) Symmetry Breaking of the HLS Model 

The HLS Lagrangian we are interested in is the lowest order expansion of Eq. (|72l ) sup- 
plemented with Eq. (1831 ). However, in order to use it with most real data, one cannot avoid 
defining an appropriate symmetry breaking mechanism. Several breaking schemes have been 
proposed E S [51 [H as there is no unique way to implement such a mechanism in the HLS 
model. We will prefer the method proposed in [[6l which looks to be the simplest that automat- 
ically fulfills the hermiticity requirement. This symmetry breaking scheme turns out to modify 
the non-anomalous Lagrangian terms in Eq. (1721 ) to : 

< 

where the SU(3) symmetry breaking 

Xa = Diag(l, 1, ./z^) 

(85) 

Xv = Diag(l, 1, y/z^) . 

As the parameter za = [fx/ fn]'^ is fixed here by kaon decay data, it can hardly be con- 
sidered as a truly free parameter, even if one allows it to vary within errors ll54l : za = 
1.495 ± 0.031. As shown in the Lagrangian pieces given in the main text, the second breaking 
parameter, zy, allows one to shift the (p meson mass away from those of the p and u mesons ; 
practically we have more freedom in varying it. The full SU(3) broken HLS Lagrangian pro- 
duced by this mechanism (without interaction terms) has been given in [16]. Using this 

^^In the Appendix 4, the two breaking parameters l^y and ix (denoted zt in the present paper) have been found 
in this reference to fulfill iiYl"^ = 1. 



-^Tr[(L - R)XAr 
-^Ti[L + R)Xv]' 



matrices Xa and Xy can be written 



(84) 
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mechanism, however, the pseudoscalar kinetic energy term of the HLS Lagrangian is no longer 
canonical and a renormalization of the pseuscalar fields is required : 

P' = X]/^PX]/^ (86) 

where P and P' stand respectively for the bare and renormalized pseudoscalar field matrices. 
With this redefinition, the kinetic energy term of the SU(3) broken Lagrangian is once again 
canonical. However, the coupling constants to kaons have to be changed correspondingly by 
introducing renormalized fields. 



With this symmetry breaking mechanism, the realm of practical relevance for the HLS 
model extends to pions and kaons, as far as pseudoscalar mesons are concerned. In order to 
bring r] and r]' mesons into the game, one needs first to define them in terms of the 773 and % 
fields, second to extend the breaking scheme from SU(3) to U(3). 
We use the one angle traditional mixing expression : 



" 7] 













cosOp —sinOp 



sin 6 1 



cos Up 



(87) 



It has been shown flV\ that the one-angle description was equivalent at the first two leading 
orders to the two-angle, two-decay constant description in favor since the Extended Chiral 
Perturbation Theory (EChPT) [fBlffllTSll. 

Now the question is how Nonet Symmetry Breaking (NSB) can be incorporated within the 
(SU(3) broken) HLS Lagrangian already defined. This can be done by means of determinant 
terms fil2il which break the f/A(l) symmetry : 



2 £2 £2 

>^^HLS + In det U ■ In det + In det d^U ■ In det ^^'U^ 



(88) 



where U is defined by Eqs.(|69l) and (1701) after removal of the a field matrix. This can be 
rewritten more explicitly : 



(89) 



Therefore, in this manner, one provides both a mass to the singlet and a modification of 
the kinetic singlet term which is thus no longer canonical and, then, calls for a renormalization. 
The exact renormalization relation is given in [1 IJ, where it has also been shown that, at leading 
order, this transformation is equivalent to using the HLS Lagrangian but replacing Eq. (|86l) by : 



P' + xP' 



X'J\Ps + P,)XT 



(90) 



(with obvious notations). The Nonet Symmetry Breaking (NSB) mechanism introduces a pa- 
rameter X which can be related [[TT| with A by : 



X 



A 2 



1 



VI + A52 



A - 0.20 - 0.25, 



(91) 
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(B = {2za + l)l'izA) with a precision better than ~ 5%. 

Therefore one has only to equip the SU(3) broken HLS Lagrangian with the U(3) broken 
renormalization condition given by Eq.(l90l). Ref. [fTTIl showed that, at leading order in breaking 
parameters one recovers the ChPT expectations. 

In order to achieve this general presentation of the broken HLS model, we recall in 
a few words the breaking procedure of the anomalous Lagrangian. A priori, the transforma- 
tion to renormalized fields given by Eq. (l90l) induces a breaking mechanism into the anoma- 
lous HLS Lagrangian given by Eqs. (|80l) and (l82l) . It has been shown [8, 9] that, alone, 
this breaking scheme (as well as no breaking at all, both) implies that the coupling constant 
ratio GK*OKO'y/GK*±K±'y equals 0.5 in sharp disagreement with experimental data [54] . In- 
terestingly, the non-relativistic quark model (NRQM) allows more freedom by exhibiting a 
dependence of this ratio upon the ratio of quark magnetic moments r [98] : 

Gk-oro^ _ 1 + ^ 

In [mm, it has been shown that this effect can be obtained by mixing a symmetry breaking 
scheme proposed by Bramon, Grau and Pancheri (BGP) [(H [51 with some sort of vector field 
renormalization. In the following Appendix, we show that the guess expressed in [^ that 
a vector field renormalization should take place is justified while having to perform the SU(3) 
symmetry breaking of the Yang-Mills piece in the HLS Lagrangian. 

Numerical analysis implies that these two mechanisms (namely the BGP mechanism and 
the vector field renormalization) are highly correlateci^ with the neat result that the broken 
VVP Lagrangian in Eq. ([83l) becomes : 

C = Ct^'P^TilXTd^i^eQA, + gV,)XT^dp{eQA^ + (?K)XtP]. (93) 

with P being replaced by renormalized fields using Eq. (|90l ) above and V being understood as 
already renormalized (as a consequence of breaking the SU(3) symmetry of Cym)- 

Therefore, the Lagrangian we use in order to account for anomalous decays i^ Eq. ([93] ) 
with : 

Xt = Diag(l, 1, v^) (94) 

where is a parameter to be fitted. We should stress that this specific breaking, which al- 
lows one to recover Eq. ([92]) leaves all other couplings of physical interest (AVP and AAP) 
unchanged. One should note that, except for (conceptually unavoidable) mixing angles, the 
model we use introduces only two parameters za and zt in the anomalous sector, the former 
being essentially fixed by pure kaon physics. Taking into account that our broken anomalous 
Lagrangian aims at accounting for 14 decay modes, the parameter freedom is actually very 
limited. 

Actually, as stated in [91, imposing that the two-photon decay widths of the 77 and 77' mesons should remain 
as given by the original Wess-Zumino-Witten Lagrangian ll95]l96l provides the functional correlation first found 
numerically. 

Eq. ( [82b might have also to be broken similarly. However, existing data on box anomalies allow access only 
to the limited sector ir^ /rj/rj' — > 7r+7r^7 not affected by more breaking than reflected by Eq. ( [90b ; it has been 
shown ||85J that Eq. ( [82] | as it already stands, suffices for satisfactorily accounting for the data. There is therefore 
no need to go beyond for this term. 
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D The Yang-Mills Term Of The HLS Lagrangian 



The HLS Lagrangian defined in the body of the text should be understood supplemented 
with the Yang-Mills piece associated with the vector meson fields [2J. Defining the vector field 
matrix as given in Eq. (l75l) (see Appendix A), this writes [2J : 



YM 



(95) 



The square of the abelian part provides the usual kinetic energy term of the vector meson 
fields and is certainly canonical while working with ideal fields. When performing the field 
transformation defined by Eqs. (fT6l) and (fTTl ). this piece remains canonical in terms of physical 
fields, as a trivial consequence of the orthogonality property of the matrix R : R(s)R{s) = 1. 

The mixed abelian-non-abelian term contributes to all vector meson self-energies. Ad- 
ditionally, it also provides transition amplitudes among the ideal p°, cu and (j) fields. One can 
check that these additional contributions to the transition amplitudes only involve K*K* loops 
and contribute in the following way to right-hand sides of Eqs. ^ : 



n 



■ ■ ■ + 9(I,k''K'9ujK''K* 

9pK*K*9ujK*K* 
np0(s) = ■■■ + gpK'K*94>K'K* 



(96) 



n 



K*±K 



still using obvious notations. Quite interestingly, one sees that the transition amplitudes n^^(s), 
^puj{s) and Yip^{s) are modified in such a way that ujcj) always receives an additional non- 
vanishing contribution, while the pu and pcj) transitions receive non-vanishing contributions 
only if isospin symmetry is broken (m/^*± ^ m^^.o). One may note the striking correspondence 
between the Yang-Mills contributions to the transition amplitudes and those already given in 
Eqs. ©. _ 

Moreover, taking into account the threshold value of K*K loops, these are certainly real 
below 1 GeV and their effects can be considered numerically absorbed in the subtraction poly- 
nomials we already use. Denoting by M the K* mass, the amputated K*K* loop fulfills : 



ImnCsl 



1 s + 3M2 (s - 4M2)3/2 



(97) 



487r M2 sV2 

above the two-i^'* threshold. This function can be algebraically derived from the two-pion 
loop (see Appendix E) and undergoes the same minimum number of subtractions as this. 
If no SU(3) breaking is implemented inside Cym, the relevant coupling constants are : 



G.. 



'K*+K* 



G 



pK*+K* 



^pK*OK'^ 



— V^G^X'^+K* 



2 

9_ 
2 



(98) 
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Finally, tadpole terms may also contribute to transition amplitudes ; they are generated by 
the non-abelian term squared. One can prove that they follow the same pattern as shown in 
Eqs. (|961) just above, namely still contributing to Iluj^{s), while contributing to np^(s) and 
Iip(f){s) only if isospin symmetry is broken. 

Another remark is of interest and concerns the flavor SU(3) symmetry breaking of the 
Yang-Mills piece. Indeed, as we already perform the flavor SU(3) symmetry breaking of the 
Ca and Cy pieces of the HLS Lagrangian (see Eqs. (f84l) and (f85l)). we should perform likewise 
with Cym- In consistency with our SU(3) symmetry breaking scheme, this should be done as : 

^YM = —-^''^[F'uuXymF'^''Xym] 

2 (99) 

Xym = Diag(l, 1, Zym) 

It is easy to check that, in order to restore the canonical form of the vector meson kinetic energy 
term, one has to renormalize the vector fields and define : 

V, = Xyl^'V^^^^Xyli' (100) 

This change of fields should be propagated to the anomalous Lagrangian. Doing this way and 
introducing the breaking procedure proposed by Bramon, Grau and Pancheri [5], one indeed 
ends up with Eq. (l93l) in Appendix O (for this purpose, see the relevant discussions in [HI HI) 
with : 

XT = Xylf . (101) 

In this renormalization, Eqs. (|981) remain valid with replacing there each G^j^^-^* by 
^4)K*K* I -^T- However, it is interesting to note that the single piece of information in our global 
data set which is really sensitive to this renormalization is the (charged and neutral) K* radia- 
tive decay. Actually, the vector field renormalization also affects the ideal ^ meson decay to the 
KK and e^e~ modes ; because of mixing among ideal fields, this also affects u; and p decays 
to e^e^ . However, one can easily convince oneself that the influence of this renormalization is 
absorbed numerically in the fitted value for the parameter zy. 

E Radiative and Leptonic Coupling Constants 

In order to simplify the main text we prefer to list here the coupling constants entering 
decay widths expression which will be treated in this paper. Most of them can be derived 
trivially from expressions already given in Appendix E in [j9l. 

E.l Radiative Decays 

Starting from the Lagrangian in Eq. (|93] ). and using the breaking procedure as defined by 
Eq. (|90l ). one can compute the coupling constants for all radiative decays of relevance. Let us 
define : 
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G 



3eg 



3eg 



K 



fn 



1 2 



K 



— , Sp = ep-eo (tan^o = 1/V2). 
za 

(102) 



Some VP'-y coupling constants are not affected by Isospin Symmetry breaking : 



G p±TY^-y 


I- 








G' 




1 

Zt 


GK'^K^-y = 


G' 


/i?(2- 


1 

Zt 



(103) 



The coupling constants are : 



G pjt^o^ 



Gprn"! 



Gpir]'f 



5- 



V2{1 — x) cos Sp — (2x + 1) sin Sp 
V2{1 - a;) sin 6p + {2x + 1) cos 6p 



(104) 



The other single photon radiative modes provide the following coupling constants : 

G 



( r 

GipjTiO^ 



G. 



Gujiri''y 
Ga 



-G v^(l — x) cos 6p — (1 + 2x) sin5p 
9 '- 

-G [(1 + 2x) cos 6p + V 2(1 — x) sin 6p 
G [Z{2 + x) cos 5p - V2Z(1 — a;) sin5p 
G [V2Z{1 — x) cos 5p + Z(2 + x) sin 5p 



(105) 



In order to go from ideal field couplings to physical vector field couplings, one has to use 
linear combinations of the couplings in Eqs. (|104H105l) weighted by elements of the transfor- 
mation matrix R{s) given in the body of the paper. 



E.2 P77 and F - 7 Modes 

The 2-photon decay modes are not affected by Isospin Symmetry breaking in the vector 
sector and keep their usual form within the HLS model LSI HI [HI : 



G 



a. 



em 



nil 



Gl 'j^i) 



cos Up — V 2 X smop 



3 

5 - 2Z 



sin 6*0 + V2 



3 

5 + Z 



X cos 9p 



(106) 
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As stated in the text, we actually replace this last coupling by using the world average value for 
as given in the RPP [[54l. 

Finally, in the non-anomalous sector, the leptonic decay widths of vector mesons depend 
on the HLS V — 7 couplings. For the ideal combinations, we have : 



fpn 



pii 



3 

fpn 



(107) 



P/7 



It was shown in that the pseudoscalar mixing angle is not a free parameter, but is related 
with the SU(3) breaking parameter Z{= 1/ za) and the Nonet Symmetry breaking parameter x 
by: 

tan^P = ^^^^x (108) 
with a very good accuracy. This relation is used in our fits as a constraint. 



E.3 Partial widths 

We list for completeness in this Section the expressions for the partial widths in terms of 
the coupling constants for the various cases which are examined in the text. 
The two-photon partial widths are : 



nil 



647r' 

The leptonic partial widths are : 

T{V e+e- 

The radiative widths are : 



Ana' 



3-I/V7I 



1 



967r 



3mi 



rriy — mp 
mv 



\GvPf\'^ 



where V is either of p°, cu, (p and P is either of 7r°, r], r]', and 



r(p ^ 1/7) 



327r 



mp — mr 



mp 



1 3 



\G 



VPj\ 



The decay width for a vector meson decaying to V + P is : 



(109) 



(110) 



(111) 



(112) 



T{V' VP) 



9671 



rjiy, — {my + mpYWrny, — (my — mp)^] 
my/ 



\Gv' 



V'VP\ 



(113) 



Finally, the partial width for a vector meson decaying into two pseudoscalar mesons of 
equal masses is : 

r(y ^ PP) = _LKi:ii!!if£!|G 



487r 



m 



VPP\ 



(114) 



V 
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F The Loop Functions 

The loop functions can be written quite generally as : 

n{s)=fis)Kis) + Pis) 



(115) 



where /(s) is a polynomial Q{s) divided by some power of s. The degree of the polynomial 
P{s) is fixed always at second degree and we require P(0) = 0. we have : 



' lmK{s) 

K{s) = Co + CiS + 023^ H / 



, (-5 > 

\mK{z) 



(116) 



dz 



ISO z^{z ~ s + ie) 

where sq is the (direct) threshold mass squared, while Sc is the (crossed) threshold mass 
squared. 



F.l The PP loop 

In the case of equal masses, Sc = and we have BUl : 



and the solution is 



s < : 
< s < So 
s> So : 
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G{s) 
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487r 
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2 (-s)V2 "^(so-s)i/2 + (-s) 



In 



■So 



- i-s) 



1/2 



1/2 



[So 



^3/2 
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2 



arctan , 



.1/2 



In- 



V 1^0 - s, 
s'/^ -(s- 



so 



,1/2 



il/2 + 



So) 



1/2 



tTT [S — So 



a/2 



(117) 



(118) 



The behavior of n(s) near s = is simply 0{s), and Q{s) behaves like 0{s) near the origin. 
This result coincides with the one of [|52l O. By performing more subtractions, one could 
choose to fix externally the actual s" behavior of the loop near the origin. 



F.2 The PP' Loop 

In this case, we have : 

n(s) = ^vpp - ^o)(s - sj ^ ^^^^^ 
487r 
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Let us define 

S < Sr 



■Sc < s < So : ^pis) = — (sq — — Sc)"^^^ arctan w ^ — 



s > So : 



-1 



ip{s) = — (s - So) / (s - s 



In 



1/2 



(s - sjV^ - (s - So) 

(s-S,)l/2+(s_So)V2 



(120) 



The solution for K is obtained by subtracting a polynomial in such a way that the behaviour of 
Eq. (11191 ) is not singular at origin : 



K{s) = (p{s) - [cq + CiS + C2S' 



(121) 



with (so = ml and Sc = ml) : 

Co = 



momc , mo-mc 
m ■ 



TT 



mo + mc 



1 mo-mc (mo - mj^ 
— m 

2tt mo + mc m^mc 

1 



(122) 



C2 



57r 



(mQ — mIY J "^0 ~ mc ^ ^m^ + m 



m^ml 



mo + mc 



m^ml 



The exact behaviour for 11 (s) at origin is then detemined by the choice of P(s). 



F.3 The PP' Loop In the Complex s-Plane 

The expressions in the two Subsections above give the value of the loop functions for any 
real values of s. It is interesting to know how these functions extend into the complex s-plane, 
i.e. for complex values of s. It is actually in this manner that the expressions above have been 
derived. One can check that : 



^{z) = --{z-Scf'\so-zf'Hn 



vr 



^So-zf'^+l{z-ScY'^ 
{So-zY/^-l{z-ScYI^ 



(123) 



is - up to a polynomial with real coefficients - the (single) analytic function of z real for 
Sc < z < So, having as imaginary part for real s > so, ImK(s), given by Eq. (11161) . The most 
general solution to Eq. (|116|) . is then written: 



K{z) = ^{z) + Pn{z) 



(124) 



with a polynomial, Pn{z), with real coefficients chosen in such a way that the behaviour at z = 
is the required one. The other coefficients have to be fixed by other external (renormalization) 
conditions. Eq. (11231) gives the loop function on the so-called physical sheet of the Riemann 
surface. The expression for (f{z) on the unphysical sheet close to the physical region s > so is 
obtained by a winding of 2n radians around the threshold (branch-point) s = so. 
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F.4 The Leptonic Loop 



In order to compute the photon vacuum polarization, one needs to have at one's disposal 
the analytic expression of the i^£^ loop. More precisely, one needs the ratio of this loop 
(lepton contribution to the photon self-energy) divided by s, the off-shell photon invariant 
mass. This can easily be derived from the function 11° (s) giver0 in [|99l or computed from 
detailed information from HI 0011 : 
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20 4 4(l-z)(l + 2z) 



where 
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2u \nu 
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(125) 



(126) 



and mi is the lepton mass and a is the fine structure constant. 

However, for explicit computation in a minimization code, on needs the explicit expression 
along the real axis. This is : 
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< -2 < 1 
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This function will be summed up with the (parametrized) hadronic vacuum polarization 
provided to us by M. Davier [|68l and H.Burkhardt ||69ll for, respectively the region above 
and below the 2-pion threshold. The function Ii{z) just defined is analytic and vanishes at 
2 = 0. The term of order can be derived from the function H*^^) given in [99] but is difficult 
to handle in fitting procedures. Other expressions for the function we use can be found in 
[|62lll01L Finally, even if possible in principle, we do not have the freedom of subtracting more 
the function 11(2;) as conditions at the Z boson mass for the full photon vacuum polarization 
seem to fix it to be zero If65l |66l |67l . 
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Figure 1: Schematic representation of the Feynman diagrams contributing to the pion form 
factor. Left plots (referred to as ai and 3.2) sketch the case of the pion form factor in e+e~ 
annihilations, while right plots (referred to as bi and b2) figure out the r decay. The upper 
plots show the non-resonant HLS specific diagrams, the lower plots describe the resonance 
contributions. The shaded blobs represent the photon vacuum polarization, the 'yV and WV 
transitions are dressed by s-dependent terms. 
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Figure 2: Data and fits for the pion form factor in e~^e~ timelike data. Left figure gives the 
fit cross section with, superimposed, the data from the Olya and CMD Collaborations [75J, 
the data set [76] from the DM1 Collaboration and the first (corrected) data from CMD2 [2T|. 
Right figure shows the form factor curve with superimposed all data sets collected recently at 
Novosibirsk ||27ll28ll29ll30ll . The region is commented upon in the body of the text. 




Figure 3: Data and fits for the pion form factor in r decay. Left figure shows the case for 
ALEPH data W, right figure shows the case for CLEO data fSSl]. 
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Figure 4: Left Figure shows the fit in the spacelike region close to s = together with the data 
from NA7 f59\ and Fermilab [80J ; the fit scale factors (1 .008 and 1 .006 respectively) have been 
applied. Rightside figure shows the prediction for the Pu phase shift with the Cern-Munich 
data [82] and the data from [|83l superimposed. 
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Figure 5: Residual distribution for all the e'^e" new timelike data over the whole invariant mass 
interval. The inset magnifies the p peak invariant mass region. 
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Figure 6: Distribution of fit residuals for the r data. The upper two plots give the residuals 
for ALEPH and CLEO data within the model presented. The lower plots shows the r data 
residuals when fitting without a p° — mass shift (i.e. Sm^ = is required) ; one should note 
the vanishing of the structure around the p peak in ALEPH data produced by the p° — mass 
shift. One should also note that the CLEO residual is not modified. 
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Figure 7: Matrix elements producing the neutral vector meson mixing. The functions shown 
are those given in Eqs. (|64l ) with their name recalled in each Figure. The upper part of each 
plot gives the real part of the function, the lower part, its imaginary part. 
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Figure 8: Isospin symmetry breaking effects following from the — uj — (f) mixing scheme. 
The upper plot shows the difference between and |F^(s)p normalized to the latter. 

The lower plot instead shows the p part of |F^(s) |^. 
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